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C D ■ The emergence of the first sources of hght at redshifts of 2 ~ 10 — 30 signaled the 

' transition from the simple initial state of the Universe to one of increasing complexity. 

We review recent progress in our understanding of the formation of the first stars and 
galaxies, starting with cosmological initial conditions, primordial gas cooling, and sub- 
sequent collapse and fragmentation. We emphasize the important open question of how 
the pristine gas was enriched with heavy chemical elements in the wake of the first su- 
pernovae. We conclude by discussing how the chemical abundance patterns conceivably 
allow us to probe the properties of the first stars, and allow us to test models of early 
metal enrichment. 
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I. INTRODUCTION 



One of the key goals in modern cosmology is to understand the formation of the first generations of stars and the 
assembly process of the first galaxies. With the advent of the first stars - referred to historically as Population III 
(Pop III) - the Universe was rapidly transformed into an incre asingly complex sys t em, due to the energy and heavy 
element input from stellar sources and accretin g black holes ([Barkana and Loehl . 1200 ll : H romm and Larsonl . 120041 : 
ICiardi and Ferrara . 2005 : Miralda-Escude . 2003f) . Anisotropics in the cosmic microwave background (CMB) allow us 
to probe the state of the Universe 370, 000 years after the Big Bang and provide us with some details of early structure 
formation. With the best available ground- and space-based telescopes, we can probe cosmic history all the way from 
the present-day Universe to roughly a billion years after the Big Bang. In between lies the remaining frontier, and 
the first stars and galaxies are the signposts of this early, formative epoch. 
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To simulate the build-up of the first stellar systems, we have to address the feedback from the very first stars on 
the surrounding intergalactic medium (IGM) , and the formation of the second generation of stars out of material that 
was enriched by the the first stellar generation. There are a number of reasons why addressing the feedback from the 
first stars and understanding second-generation star formation is crucial: 

• Over the past fifty years, there has been extraordinary progress in our understanding of the origin and evolution 
of the chemical elements over cosmic time. But there remains considerable uncertainty on the physical processes 
that seeded many elements in the first stars. 



The initial burst of Pop III star formation may have been r ather brief due to the st r ong negative feedback 
effects that likely acted to self-limit this formation mode ( Greif and Bromnil . 2006 : lYoshida et al . 2004i ). 
Second-generation star formation, therefore, may well have been cosmologically dominant compared to Pop III 
stars. A subset of the second-generation stars with masses % 0.8 M.q must have survived to the present day. 
These stars provide an indirect window into the Po p III era, as they refl e ct the enrichment f rom a single, or 
at niost a s mall multiple of supernova (SN) events (Beers and Christliebl. 2005 : Frebel al. . 2007 : Karlssonl 
20051 l2006t iKarlsson and Gustafsso^ . 120051: iKarlsson et aLl . l2008l : ISalvadori et a^.l . 120071: iTumlinsonl . |2006[) . Bv 



scrutinizing their chemical abundance patterns, we can extract empirical constraints which we apply to super- 
computer simulations, which in turn allow us to derive theoretical abundance yields to be compared with the 
data. 



• The first steps in the hierarchical build-up of structure provide us with a simplified laboratory for studying 
galaxy formation. Did the rcionization epoch influence the subsequent evolution of galaxies? What was the 
role of galactic winds in the first galaxies? Did star formation occur in bursts, or in a steady, self-regulated 
mode? How were the first nuclear black holes seeded? We can probe these formative times by reconstructing 
the conditions in the first galaxies from the chemical signaturefo in the most ancient stars. 

Existing and planned observatories, such as the Hubble Space Telescope (HST), the 8-lOm class telescopes and 
the James Webb Space Telescope (JWST), planned for launch in 2018, yield data on stars and quasars less than 
a billion years after the Big Bang. The ongoing Swift gamma-r ay burst (GRB) m i ssion provides us with a window 
into massive star formation at the highest observable redshifts (jBromm and Loebl . |2002| . 120061 : iLamb and Reichart . 
20001) . Measurements of the near-IR cosmic background radiation, both in terms of the spectral energy distribution 



and the angular fluctuations provide additional constraints on the overall energy production due to the first stars 
(iDwek et al . 2005t Fernandez and KomatsuL 2006; Kashlinskv et al. . 2005t Magliocchetti et al . 2003t Santos et al. 



20021) . Focusing on the nearby Universe, the Gaia missioro will provide six-dimensional phase-space information on no 



less than one billion stars in the Milky Way and projects like HERME ^, APOGEE0, LAMOStI, and the Southern 
Sky SurvejH, will measure detailed chemical abundances for millions of Galactic stars, out of which some fraction will 
likely probe the very first metal enrichment phase. Understanding the formation of the first stars and galaxies is thus 
of great interest to observational studies conducted both at high redshifts and in our local Galactic neighborhood. 

In this review, we will focus on the interplay between primordial star formation and pre-galactic me tal enrichment, 
since it provides one of the main feedback mechanisms that shaped the early I GM dMadau et al. . 2001 ). The chemical 
feedback from the first SNe had a far-reaching impact on early cosmic history (jCiardi and Ferrarall2005l ). Our present 
understanding is that the character of star formation changed from the early, high-mass dominated Pop HI mode to 
the more normal, lower-mass Pop II mode once a critical level of enrichment had been reached, the so-called critical 
metallicity, ZcHt ^ 10""* Zq (Omukai, 2000; Bromm et al, 2001; Bromm and Loeb, 2003; Schneider et al., 2003). It 
is then crucially important to understand the topology of early metal enrichment, and when a certain region in the 
Universe becomes supercritical. In general, metals produced by Pop III SNe will initially be dispersed into the IGM, 
and a fraction of them will later be re-incorporated into more massive systems during bottom-up structure formation. 
The metal-enriched protogalaxies then drive further heavy element synthesis and ejection into the intergalactic medium 



^ In this review, we make a clear distinction between a chemical fingerprint and a chemical signature (see Sec. lIVIl . The term chemical 
fingerprint is used when the elemental abundances provide direct evidence of a specific reaction mechanism: the r-process or triple-a 
process, for example. This review is concerned with identifying the chemical signatures of the first stars in the surface abundances of 
the oldest stellar populations. A chemical signature is inherently more complex because the elemental yields expelled from a dying star 
are likely to depend on more than one physical process. The signature then reflects multiple parameters, such as stellar mass, rotation, 
explosion energy, and the amount of fallback onto the remnant. 

^ http://www.rssd.esa.int/gaia 

^ http://www.aao.gov.au/hermes 

* http:/ /www. sdss3.org/surveys/apogee.php 

^ http://www.lamost.org 

® http:/ /www. mso.anu.edu.au/skymapper 
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on a more massive scale. They are also candidate drivers of the early hydrogen reionization that has been inferred 
from the relatively large optical depth in the CMB polarization measured with the Wilkinson Microwave Anisotropy 
Probe (WMAP). This second wave of metal injection pollutes a larger fraction of the Universe, and sets the stage 
for pervasive gravitational fragmentation in those dark halos that have avoided complete photoevaporation in the 
emerging cosmic ionizing background. 

We now introduce our working definition of what we mean by a "first star" . A Pop III star contains no heavy 
elements, such that its metallicity is Z — 0. A Pop II star, on the other hand, has Z > Zcrit- In principle, there could 
have been stars with < Z < Zcrit, but none of them, by the definition of critical metallicity, would have survived 
until today. More importantly, such hypothetical stars were extremely rare, if current cosmological simulations are on 
the right track. Specifically, those simulations are indicating that even a single Pop III SN can enrich its surroundings 
to values above ZcrU, or that gas with Z < Zcrit resides in low-density regions of the IGM, where star formation does 
not occur. Once we learn more about the topology of early metal-mixing, we may wish to revisit the definition for 
Pop III, but for now, our terminology should work well. For the record, no metal-free star has been observed in any 
of the major surveys that seek to identify the most metal-poor stars at the present epoch, i.e. a spectrum exhibiting 
H absorption lines with no trace metals (T. Beers, personal communication, 2012). The lack of metal- free stars fits 
within a developing paradigm in which the first stars were very massive and short-lived (see Sec. |ll]), although if such 
an object is eventually found, there will be no shortage of theories to support its existence! 

A prime objective for current astrophysics is to predict the properties of the first galaxies with the goal of detecting 
them with the JWST or an Extremely Large Telescope (ELT). What is the required minimum stellar and, by extension, 
virial mass in order for them to be detectable? This will determine the number of sources accessible through deep field 
campaigns ( Pawlik et aLl . [201l[ ). In turn, it will depend on the stellar population mix within the first galaxies. Is the 
fine-grain mixing of the metals that fall into the center of the galactic potential well incomplete, so that Pop III stars 
would form simultaneously with Pop I/II stars? Or is the turbulent mixing so efficient as to wipe out all remaining 
pockets of primordial gas? Thus, arguably, the most important ingredient that defines the physical state of the first 
galaxies is their metal content. This is clearly a very complex problem. Recently, it has become possible to attack it 
with numerical simulations, utilizing the latest breakthroughs in algorithm development and supercomputing power. 

Here, we raise a number of issues that are highly controversial in the general literature, and appear in need of 
clarification. Because of recent developments in galactic dynamics, we argue that it is not obvious where astronomers 
should look for the most ancient stars. Pub hshed numerical simulations are already unclear on whether the o ldest stars 
are solely the preserve of the inner bulge ( Bland-Hawthorn a.nd Peebled . l200fil : IWhite and Springe! 1200(11) or spread 

Scannapieco et alT 20061) . Given that galaxies appear to grow inside-out 



over the entire Galaxv (iBrook et al . 2007 



( Cooper et ai . 2010t Zolotov et all 2009 ) , it is often a s sumed that the majority of these stars must be confined to 
the inner bulge, consistent with its age (jZoccali et all 120061 ). But regardless of the galaxy formation model, the 



long-term dynamical evolution of the galaxy must also be considered, and this has important consequences for our 
understanding of the early Milky Way (Sec. IV. B|) . 

The uncertainty in the birth site, and therefore the birth process, leads us to question another tenet of stellar 
astrophysics. Does it necessarily follow that targeting the most metal-poor stars is an efficient route to learning about 
the yields of the first stars? This needs some clarification. A low valu£0 of [Fe/H] is no guarantee that a star is 
ancient since it may reflect environmental conditions (e.g. low star formation efficiency, shallow potential well). But 
radioactivity in som e extremely metal poor stars provides clear independent evidence that many are indeed ancient 
(ISneden et aLl . l2008l) . 

Conversely, a high value of [Fe/H] does not indicate that a star is of young or intermediate age. Super-solar 
abundances are detected in sources out to the highest detectable redshifts because the dynamical times are very s hort 
in the cores of galaxies ( Freeman and Bland-Hawthorn , 2002 ; Hamann and Ferland , 1999t ISavaglio et al. , 2012 ) . It 
is now well established that the bulge, the halo and all dwarf galaxies show a spread in [Fe/H] and comprise stellar 
populations that are 10 Gyr or older, equivalent to a formation redshift prior to z ^ 2. Globular clusters, which are 
some of the most anc ient systems with well dete rmined ages, exhibit a bimodal distribution in metallicity in the range 
—2.5 < [Fe/H] < ( Brodie and Stradeii 120061 ). These are unlikely to retain a clear chemical signature of the first 
generation of stars ( Larsen et all 20121) . 

Given the large errors in the ages of old stars, the relative fractions of these stars that formed before, during or 
after the reionization epoch is an open question. In other words, there are no obvious stellar age-metallicity relations 
that we can appeal to at the present time to help us unravel the sequence of events in the early universe. We believe 



[A/B] = logig(nA/nB)* — logio('*A/'T'B)0, where nx is the number density of element X. 
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that this complex situation can only be sorted out with a far greater understanding of the origin and evolution of the 
chemical elements. 

Our working definition of what constitutes a "first star" needs further clarification regarding the timing of star 
formation. A relatively small number of widely dispersed stars (and conceivably star clusters) triggered the reionization 
epoch in the redshift interval z ~ 15 — 20. Evidently, these stars belong to the first stellar population, i.e. to Pop III. 
But we can extend the temporal definition of a "first star" to include those that were the first to enrich their immediate 
environment regardless of epoch. Most of the gas confined by collapsing dark matter will not have experienced any 
form of chemical enrichment. The stars that formed in these regions during or after the reionization epoch are the first 
to enrich the local gas. This stellar population is identified by location rather than by epoch but since it specifically is 
formed out of virtually metal-free gas, it too is ascribed to Pop III. As we shall see, however, due to the all-pervasive 
ionizing field, the two sub-populations are thought to have been physically distinct. 

It is not known at what stage in cosmic time the last vestiges of pristine gas finally succumbed to stellar enrichment, 
but it is not inconceivable that some Pop III stars formed a billion years after the reionization epoch. This raises 
the tantal izing prospect that we can identify s uch re gions in direct observations of the intermediate and high-redshift 



Universe ( Johnson . 2010t Scannapieco et al. . 2003 ). We anticipate that Pop III stars will have distinct chemical 



signatures that can be identified in large stellar surveys. 

We direct t he re ader to excellent reviews on early star formation from which we have drawn inspiration. 
iBromm et ali (120091 ) discuss the state-of-the-art numeri cal simulations of primordial star formation and the as- 
sembly of the first galaxies. iBeers and Christliebl (l2005h disc uss the search for, and detailed observations of the 
most-metal poor stars in the Galaxy, whi le iTolstov et air(|2009f l focus on the kinematics and chemical abundances of 
stars in dwarf galaxies, and iHelmil ( 2008 ) describes the formation of the Galactic stellar halo and its connection to 
dwarf galaxies. In this review, we focus on the chemical signatures of the first generations of SNe. These signatures 
provide us with an important empirical prob e to test our t heories of hi gh-redshift star formation with o bservations of 
metal-poor stars in our cosmic neighborhood ( Frebel et aL ,i 2009; Karlsson et al . 2008 : Tumlinsonl . 2006[ ). an approach 
often called Stellar Archaeology or Near-Field Cosmology ( Bland-Hawthorn and FreemanL |2000[ ) . 



II. FORMATION OF THE FIRST STARS 



Th e first stars in the Universe likely formed roughly 200 Myr after the Big Bang (jAbel et all l2002t iBromm et al. 



19991 . |2002| ). when the primordial gas was first able to cool and collapse into dark matter minihalos with masses of 
the order of 10^ Mq (see Fig. [T]). These stars are believed to have been predominantly very massive, with masses of 
~ 30 — 1OOA^0, owing to the limited cooling ability of primordial gas, which in minihalos could only cool through 
the radiation from H2 molecules. While the initial conditio ns for the forma t ion of these stars are, in principle, known 



from precision measurements of cosmological parameters (jKomatsu et al\ . 120111 ). Pop III star formation may have 
occurred in different environments which may have allowed for different modes of star formation. Indeed, it has 
become evident that Pop III star formation may actually consist of two distinct modes: one where the primordial gas 
collapses into a dark matter (DM) minih alo, and one where the metal- free gas becomes significantly ionized prior to the 
onset of gravitational runaway c ollapse ( Johnson and Bromm ^ 2006 ). This latter mode of primordial star formation 
was originally termed 'Pop II. 5' (jGreif and BrommT 2006H Johnson and BrommL [2OO6 : iMackev et all 20031 ). To more 
clearly indicate that both modes pertain to metal-free star formation, we here follow the new classification scheme 
( Johnson et al I I2OO8I : lMcKe7 and Tanl |2008|) . in which the minihalo Pop III mode is termed Pop III.l, whereas the 
second mode (formerly 'Pop II. 5') is now called Pop III. 2. 



A. Population III: The first mode 



The formation of the very first stars. Pop III.l in the new terminology, can be understood with two basic ingredients: 
ACDM structure formation and the atomic and molecular physics of the primordial hydrogen and helium. Here, 
ACDM refers to a Universe composed of cold dark matter (CDM), but dominated by dark energy, possibly Einstein's 
cosmological constant (A). It has been recognized for some time th at within variants of CDM, minihalos provide the 



first star formi n g site s, where cooling relies on molecular hydrogen ( Couchman and Reed . 1986 : Haiman et al . 119961: 



Tegmark et al . 1997t ). The properties of the dark matter must be drastically modified before one sees significant 



deviations fr om this robust basic prediction, such as assuming warm dark matter (WDM) models, or s elf- annihilating 
dark matter ( Gao and Theunsl . [2007HRipamonti et aLl . l2007HSpolvar "oZI . l2008HYoshida et allkoO^ . The question 
then is: What kind of stars emerge during the H2-facilitated collapse of the pure H-He gas into the minihalo DM 
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FIG. 1 The first stars form in minilialos. Sliown is tlie projected density on progressively smaller scales, reaching from the 
entire simulation box (top left) down to the center of the newly virialized minihalo on scales of 10 pc (bottom left). In the 
bottom two panels, the asterisk denotes the location of the first sink formed. The numerical technique employed here uses 
sink particles as proxies for growing protostars. Notice how the morphology approaches an increasingly smooth distribution, 
corresponding to the characteristic loitering state just prior to gravitational runaway collapse. From Stacy et al. (2010). 



potential wells? This problem, although still beyond our observational horizon, is much simpler than the corresponding 
star formation process in the local, well observed Universe where the environment of giant molecular clouds is extremely 
complex. Indeed, prior to the formation of the first stars, the Universe had no heavy elements, and therefore no dust 
to complicate the physics of cooling and opacity. It was also likely that magnetic fields did not yet play a dynamically 
significant role (but see Schleicher et ai, 2009, 2010). Finally, the early post-recombination Universe was devoid of 
external ionizing radiation fields, and strong drivers of turbulence, such that the collapse into the minihalos proceeded 
in a rather quiescent fashion. The formation of the first stars inside minihalos, where we know the initial conditions 
as given by ACDM cosmology, and where we have a complete understanding of the relevant physical processes, thus 
provide s us with a well-posed problem, amenable to rigorous numerical studies. Since the late 1990s, a number of 

bromm e^"aIl . ll999ir2002i : lNakamura and Umemural - boOlHO'Shea and Normanl . 12007 . 



provide s us with a well-posed 
group s (I Abel et aLl.l20od.l2002t 



20081 : lYoshida et all 120061 ) have simulated the formation of the first stars with sophisticated numerical algorithms, 
utilizing either smoothed-particle hydrodynamics (SPH), or adaptive-mesh refinement (AMR) techniques. These 
calculations have converged on a number of main results, leading to the current 'standard model' of first star formation, 
although important open questions remain (see below). 

The most important result, where there is general agreement, is that Pop III.l stars were predominantly massive. 
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To first order, this can be understood as the consequence of a large Jeans mass in gas that cools only via H2. Prior to 
under going runaway collap s e, the primordial gas settles into what is sometimes termed a quasi-hydrostatic 'loitering' 
state (jBromm and LarsonL 120041 ) . This state is characterized by typical values for the temperature and number 
density: T^har — 200 K and richar — 10'*cm~^. Recent numerical simulations indicate that Pop III.l stars formed with 
characteristic masses of > 30 Mq. It is likely that these stars formed with a range of masses, described by the 
so-called initial mass function (IMF). The IMF gives the number of stars formed per unit mass, where present-day 
star formation is often described with a power-law: dN/dm^ cx 'm~^ , or a sequence of such power- laws (Kroupa 
2001). For Pop III, we do not know the complete functional form of the underlying IMF with any certainty. There 
is, however, a simpler, first-order way to characterize the outcome of the star formation process by focusing on the 
mass-average: M^, oc / rrii^dN. This mass is what we mean by "characteristic mass" as it describes the typical result 
of star formation, with the understanding that some stars form with lower, and some with higher masses. It has, 
however, not yet been possible to self-consistently simulate the assembly of an entire Pop III sta r, starting from 



realistic cosmological initial conditions. Recently, such ab initio calculations (jYoshida et all |2008a|) have traced the 



evolution up to the point where a small protostellar core has formed at the center of a minihalo. This initial hydrostatic 
core has a mass, m^, ^ 10^^ ■Mq, very similar to present-day, Pop I, protostellar seeds. The s ubsequent growth of 
the p rotostar through accretion, however, is believed to proceed in a markedly different way (jBromm and Larson , 
2OOJ). In the early Universe, protostellar accretion rates are believed to have been much larger, due to the higher 



temperatures in the star forming clouds, which in turn is a consequence of the limited ability of the primordial gas to 
cool below the ^ 100 K accessible to H2-cooling. The higher accretion rates, together with the absence of dust grains 
and the correspondingly reduced radiation pressure that could in principle shut off the accretion, conspire to yield 
heavier final stars in the early Universe. Estimates for the masses thus built up arc somewhat uncertain, but a rough 
range is ^ 30 — 100 M.q- It is possible that this range may extend to lower, as well as higher values, depending 
on the details of the still very uncertain primordial IMF. The current frontier in numerical simulations attempts 
to carry out fully self-consistent radiation-hydrodynamical calculations of the Pop III accretion process, taking into 
account effects of negat ive protostellar feedback and of centrifugal support in the accretion disk around the first stars 
(|McKee and Tanl . I2OO8I ) . First attempts have now been made (Hosokawa et al. 2011; Stacy et al. 2012), suggesting 
the range of masses as given above. 

Although such fully self-consistent simulations are not yet available, an import ant new development is that Pop III. l 
stars may have typically formed as members of a binary or small multiple system ( Stacv et al. . 20 id : Turk et al. . 2009f l. 
In simulations with sink particles, where the evolution can be followed beyond the initial collapse of th e first high 



density peak, one sees the emergence of a compact disk around the first protostar ( Clark et al. . 2011bl ). This disk 



is gravitationally unstable and fragments into a dominant binary, possibly with a few more lower-mass companions 
(see Fig.[2|). These simulations have not yet been able to reach the asymptotic end state, where all available mass is 
either accreted or permanently expelled. The reason is again that radiative feedback processes in partially optically 
thick mat erial cannot be neglec t ed once a protosta r has grown to IOA^q, when full radiation- hydro simulations are 
required (jHosokawa et al. I, I2OIII : IStacv et al\ . l2012f) . Currently, the binary statistics of Pop III, regarding quantities 
like mass ratios, orbital parameters, and overall binary fraction are not yet known with any certainty. Rapid progress 
is to be expected, however. 



B. Population III: The second mode 



While the very first Pop III stars (so-called Pop III.l), with m asses of th e order of 100 J^c^, forme d within 



DM minihalos in which primordial gas cools by H2 molecules alone (jAbel et all 120021 : iBromm et all l2002f ) , the HD 



molecule can play an important role in the cooling of primordial gas in situations where the gas experiences substantial 
ionization. The temperature can then drop well below 200 K. In turn, this effi cient cooling may lead to the formation 



of primordial stars with characteristic masses of the order of ~ 10 Al© ([Johnson and Bromml . |200i), so-called 
Pop III. 2 stars. In general, the formation of HD, and the concomitant cooling that it provides, is found to occur 
efficiently in primordial gas which is strongly ionized. This is largely due to the high a bundance of electrons which 
serve as catalyst for molecule formation in the early Universe ( Shapiro and Kan3 . 1987 ). 



Efficient cooling by HD can be triggered within the relic H 11 regions that surround Pop III.l stars at the end of their 



brief hfetimes, owing to the high electron fraction that p ersists in the gas as it cools and recombines ([Johnson et al. 



2007; Nagakura and Omukai . 2005t IVoshida et~al. . 2007a). The efficient formation of HD can also take place when the 



primordial gas is collisionally i onized, such as behind the shocks driven by the first SNe or in the virialization of niassiv e 
DM halos ( Greif and Bromm . 2006 : Johnson and Bromm . 2006 : Machida et al. . 2005 : Shchekinov and Vasilievl . 12006) . 



There is a critical HD fraction, necessary to allow the primordial gas to cool to the temperature floor set by the CMB 
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FIG. 2 Kinematics in the disk-like configuration in the center of a minihalo (cf. Fig. [T]). Asterisks denote the location of the 
most massive sink, crosses of the second most massive sink, and diamonds of the other sinks. Top: Density projection in the 
X — z plane after 5000 yr, shown together with the velocity field. Velocities are measured with respect to the center of mass 
of the gas distribution. Bottom: Same as above, but in the orthogonal (x — y) plane. It can be seen how an ordered, nearly 
Keplerian, velocity structure has been established within the disk. From Stacy et al. (2010). 



at high redshifts: Xhd = nun/n ^ 10^^, where nno is the number density of HD molecules and n that of all particles. 
Except for the gas collapsing into the virtually un-ionized minihalos, the fraction of HD typically increases quickly 
enough to play an important role in the cooling of the gas, allowing the formati on of Pop III.2 stars. A n interesting 
environment for Pop III.2 formation may be the so-called atomic cooling halos ( Oh and Haimanl . |2002| ). They have 
been suggested as candidates for the fi rst galaxies to form (see Bromm et al., 2009). If these systems manage to 



20081 ) , the HD cooling channel will become important in tying temperatures to the 



remain metal-free (jJohnson et al. 
CMB (see Fig. [3]). 

There may thus be a progression in the characteristic masses of the various stellar populations that form in the early 
Universe. In the wake of Pop HI. I stars formed in DM minihalos, where > 30A^q, Pop III. 2 star formation, where 
Mi, ~ 1OA^0, ensues in regions which have been previously ionized, typically associated with relic H ii regions left over 
from massive Pop III.l stars collapsi ng to black holes, while even later, when the primord ial gas is locally enriched with 
metals, Pop II stars begin to form (jBromm and Loebl . l2003t iGreif and Bromm . l2006l) . Recent simulations confirm 
this picture, as Pop III. 2 star formation ensues in relic H II regions in well under a Hubble time, while the formation 
of Pop II stars after the first SN explosions is delayed by more than a Hubble time (Grcif et al., 2007; Yoshida et 
al.. 2007a, b; but see Whalen e t al, 2008). A nother key question is related to the role of turbulence in shaping the 



primordial mass function (|Clark et al. . 2011a[ ). If sufficient turbulence were present, a broad range of fragment masses 
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a B 21, .1 m to 

FIG. 3 Cooling channels in primordial gas. Gas properties inside a minihalo (left-hand panel) and an atomic cooling halo (right- 
hand panel). Shown are the temperature, electron fraction, HD fraction and H2 fraction as a function of hydrogen number 
density, clockwise from top left to bottom left. Left-hand panel: In the minihalo, adiabatic collapse drives the temperature to 
~ 10^ K and the density to nn ~ 1 cm~^, where molecule formation sets in allowing the gas to cool to 200 K. At this stage, the 
central clump becomes gravitationally unstable and eventually forms a Pop III.l star. Right-hand panel: In the first galaxy, a 
second cooling channel has emerged due to an elevated electron fraction at the virial shock. The molecule fraction is enhanced 
in turn, thus enabling the gas to cool to the temperature of the CMB. Under these conditions. Pop III. 2 stars are predicted to 
form. The solid green lines denote the path of a representative fluid element that follows the Pop III. 2 channel. From Greif et 
al. (2008). 



could result, similar to the power-law extension towards high masses, observed in the present-day IMF. 



III. CHEMICAL FEEDBACK 

In this section, we discuss the transport and mixing of metals in the gaseous medium between the stars. The 
efficiency of this mixing is crucial to how, where, and when the transition from Pop III to Pop II star formation 
occurred. We begin by briefly reviewing the final fate of metal-free stars, stellar nucleosynthesis, and the current 
status of stellar yield calculations. 



A. Final fate of the first stars 



Unless a star loses a significant fraction of its mass in a stellar wind, the final fate of the star is predominantly 
determined by its initial mass (see Fig. |4|). Generally speaking, this holds true also for the nucleosynthesis in the 
star's interior. Primordial massive stars should generally have weak stellar winds due to lower atmospheric opacities 
but this depends critic ally on whether the outer layers essentially remai n metal-free during the lifetime of the star 
( Ekstrom et ^jikewise, mass loss driven by nuclear pulsa t ions (ISchwarzschild and Harml . I1959I) : or metal- 

free, very massive stars is probably low in most cases (jBaraffe et al. , 2001 ) . 

Table U shows, for non-rotating stellar models (however, see Ohkubo et al., 2006 for rotating 500 and 1000 A4q 
models), the predicted final fates of massive primordial stars and the remnants that are left behind, for different mass 
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FIG. 4 Initial-final mass function of non-rotating primordial stars. Notice that there are distinct mass regimes where enriched 
mate rial is expected to be ejecte d into the surroundings, separated by ranges where no or very little enrichment would result. 
From lHeeer and WooslevI (|2(302l l. 



regimes. The various ranges of initial stellar masses should be taken as indicative as they may change if properties like 



rotation and mag netic fields are shown to be significant for the evolution of metal- free stars ( Ekstrom et ai . 2008bl: 



Stacv et al\ . \20m 



At the end of their lives, stars below ~ 9 Mq commonly undergo a series of deep mixing events where He burning 
products, typically carbon, are dred ged up to the s i irface . In particular, for metal-free stars below ^ 3.5 TW©, this 
process is triggered by a mechanism (jFuiimoto et al. I. ll990l) called helium flash-driven deep mixing (He-FDDM) rather 
than, or in addition to, the third dredge-up (TDU), operating in metal- rich, asymptotic giant branch (AGB) stars. 
Primordial stars > 3.5 A^©, may experience TDU, bu t no He-FDDM, if convective oversho ot is able to inject freshly 



20011: ISiess et all 120021 ) . With carbon mixed into 



synthesized carbon into the hydrogen-rich outer layer (jChieffi et 
the a tmosphere, metal-free intermediate-mass stars should be able to develop a dust-driven wind (jMattsson et 



20081 ) in a similar fashion to their metal-rich cousins. Consequently, the entire envelope would eventually be blown 



away and leave behind a white dwarf. 

Stars in the mass range 9 < mj M.q < 40 are believed to explode as core collapse SNe. At the end of hydrostatic 
burning, stars above ^ 10 M.q contain an electron degenerate iron core. As the Fe core continues to grow by silicon 
shell burning, the electron pressure can eventually no longer counterbalance the increasing gravitational pull and 
the core collapses. An explosive instability develops when the infalling outer layers bounce off the collapsed core 
(proto-neutron star) a nd are ejected i nto the interstellar medi um (ISM). The details of the explosion mechanism is 
still a topic of debate (jBurrows et "all. l2006t Ijanfa" et a/.l . I2OOII) . Stars around 9 — 10 M.q instead form a degenerate 
0+Ne-f-Mg core which may collapse due to rapid electron captures on ^"Ne and ^''Mg, prior to the ignition of Ne 
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TABLE I Final fates and remnants of massive primordial stars for different initial masses. 



Mass range 


Final fate 


Metal 


Remnant 


{Me) 




ejection 




m< 9 


Planetary nebula^ 


Yes 


White dwarf 


9 < m < 10 


O/Ne/Mg core collapse SN'' 


Yes 


Neutron star 


10 < m < 25 


Fe core collapse SN^ 


Yes 


Neutron star 


25 < m < 40 


Weak Fe core collapse SN"* (fallback) 


Yes 


Black hole 


25 < m < 40 


Hypernova^ (fallback) 


Yes 


Black hole 


40 < m < 100 


Direct collapse'' (no SN) 


No 


Black hole 


100 < m < 140 


Pulsational pair-instability SN^ (fallback) 


No? 


Black hole 


140 < m < 260 


Pair-instability SN'' 


Yes 


No remnant 


260 < m < 10^ 


Direct collapse_^core coUapse'^'^ (no SN?) 


No? 


Black hole 


10^ < m 


Direct collapse before reaching main sequence'* 


No 


Black hole 



"'Hcgcr et aV. ('2003). 

Nomoto et al. (200f|). 

Fryer et al. (2001). 
''lOhkubo et ai. 120061). 



( Barkat et al. 1. 119741: iNomotd . I1981T) . If s o, an electron-capture, or ON eMg SN is formed. Alternatively, the star loses 
its outer layers and forms a white dwarf (jGarcia-Berro and Ibenl . 119941 ) . At the high mass end of the core collapse SN 
regime, 25 — 40 M.q, the potential energy of the stellar envelope is comparable to the kinetic energy of the explosion 
and the inne rmost layers of the star fall back onto the central neutron star which eventually collapses to a black hole 

Tigj) . 



(|Frvei . 



A weak, or faint SN is formed if the black hole is non - rotating. For rotating blac k holes, however. 



, , 

a much stronger explosion may instead be expected (iFrver and Heger , 2000t Nomoto et al}. '200 6|) in the f orm of a 
hypernova (HN). SN1998bw (jCalama ei"aZl . ll998Hrwamoto et a/.l , 119981 ) and SN20031w (|Mazzali e t all l2006f ) are two 
examples of unusually energetic SNe commonly ascribed to the HN branch. Stars with masses above AO Mq collapse 
directly to a black hole (see, however, Ekstrom et al, 2008a). 

For Pop. Ill stars with masses ^ 100 A^o, the pair instability kicks in after central carbon burning, driven by 
the creation of electron-positron pairs. Below ~ 140 Mq, this instability causes the star to pulsate violently. While 
the outer layers are lost in SN-like explosions, the star settles down to form a massive Fe core. Eventually, the star 
collapses quietly to form a black hole. On the other hand, if the star is in the mass range 1 40 < m/M.^, < 260 
it wi ll face complete disruption as a result of the pair instability, and no remnant is formed (jHeger and Wooslev , 
l2002l) . In still more massive stars, the energy released from explosive oxygen and silicon burning, caused by the pair 
instability, is instead used to photodisintegrate the nuclei in the c entral core. The explosion is halted and a black hole 
is formed (jFrver et al ]. l200l[ ). As discussed bv lFryer et al\ (|200l[ ). rotation could, however, change this scenario (see 
also Ohkubo et al., 2006) and stars above ~ 300 M.q may be able to explode as SNe, in which tremendous amounts 
of energy exceeding 10^'' erg would be released. In contrast, supermassive stars ^ 10^ -M©, if th ese exist, will collapse 
even before settling down onto the main sequence, ow ing to a general relativistic instability (|Ohkubo et al. .1. l2005 



Shapiro and Shibatal . 120021: IShibata and Shapirol . 120021 ). 



B. Nucleosynthesis in the first stars 

Whether the primordial stars will eject any metals into the ISM at the end of their lives is of utmost importance 
for the early cosmic chemical enrichment, as well as for how subsequent (Pop. II) star formation proceeds. In Sec. 
llVi we shall discuss in detail which chemical signatures are observed and what these may be telling us about the Pop. 
Ill stars and the primordial IMF. In preparation for that section, we will here briefly review what possible chemical 
fingerprints and signatures to expect from primordial stars of various masses. Table U shows in which mass regimes 
we expect the stars to enrich their surroundings in metals and in which regimes they do not. In principle, stars below 
rougly 40 A^Q will be able to eject freshly synthesized material, through SNe or winds. The same goes for stars in 
the mass range 140 < mj M.q < 260, which are completely disrupted as they explode as pair-instability SNe (PlSNe). 
Stars in the range 40 < m/M.Q ^ 100 are expected to undergo direct collapse without an explosi on, while stars in the 
pulsa tional pair-instability SN regime presumably do not eject any elements heavier than helium ( Heger and Wooslevl . 



20021 ). These two groups are therefore omitted in the detailed discussion below. Similarly, stars in the mass range 
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TABLE II Nucleosynthetic signatures of massive primordial stars. 



Mass ran, 




Expl. energy 


Ejected Fe 


Metal 


Nucleosynthetic 






(xlO^^ ere) 




GnrichmGiit 


characteristics 




9 


Wind 





^Li(?), C, N, O, 


(C-hN)/0 > 1 










Na, Mg, s-process 


e.g., [Pb/Ba] ~ 1.2 












[C, a/Fe] < 


9 < m < 


10 


^ 0.1 


- 0.002 - 0.004 


Carbon to Fe-peak^ 


[Mg/Ca] < 










r-process? 


[Ni, Zn/Fe] > 












e.g., [Ba/Eu] ~ -0.6^ 


10 < m < 


25 


~ 1 


- 0.07 


Carbon to Fe-peak 


[a/Fe] > 


25 < m < 


40 


^ 1 
< 1 


^ U.Ul 


Carbon to Fe-peak 


[C, O/Fe] > 










r-process? 


e.g., [Ba/Eu] ~ -0.6*' 












larger [Si,S/C,0] 


25 < m < 


40 


>io 


~ 0.08 -0.3 


Carbon to Fe-peak 


larger [V,Co,Cu,Zn/Fe] 












smaller [Mn,Cr/Fc] 


40 < m < 100 


No expl. 





— 


— 


100 < m < 140 


~ 1 





Only H, He(?) 














[Mo- m /Nn A 11 ^ n 


140 < m < 260 


~ 10 - 10^ 


~ 0.01 - 40 


Carbon to Fe-peak 


[Si,S/C] ~ 1 - 1.5 












[Zn/Fe] <C 0, no r-proc. 












[Mg, Si/Na, Al] > 


260 < m < 


lO'^ 


~ 10^ - 10" 


~ 5 - 20 


Carbon to Fe-peak_^ 


[C,Q/Fe] < 












[Zn/Fe] > 


10" <m 




No expl. 










" Based on a solar metallicity 8.8 Mq ID-model llWanaio et aLl. l2009l) . 

^ Deduced from observations of metal-poor stars ("Barklem et al., 2005). 

Ejection of metals may only occur if bipolar jets are generated (,Ohkubo et aZ.L 120061) . 



260 < m/ M.Q < 10^ are expected not to eject any metals, at least not in non-rotating models, as everything would be 
swallowed directly by the central black hole before an explosion takes place. Howe ver, in the case of rot ation, bipolar 
outflows may develop and significant amounts of metals could possibly be ejected ( Ohkubo et aZ.I . [2006() . 

It should be noted that current SN yield calculations are performed without knowledge of the proper explosion 
mechanism(s). In the literature, various methods such as a piston, energy injection, or enhanced neutrino opacity, are 
used to drive the explosion (see Fryer et al., 2008 for a discussion). Moreover, due to the unknown amount of angular 
momentum in the stellar core, the fraction of mass falling back onto the proto-neutron star is uncertain and, like 
the explosion energy and the electron fraction Y^, is treated as a model parameter. Consequently, yields of elements 
synthesized in the layers closest to the center of the star, e.g. the Fe-peak, are marred by significant uncertainties. 
Only when the explosio n can be treated in a self-co nsistent way can these issues be resolved. Most models also lack a 
description of rotation. iMevnet and Maedeii (|2002l ) have shown that the surface abundance of ^''N is greatly enhanced 
in metal-poor, massive stars as a result of rotationally induced mixing (see also Meynet et al., 2006; Hirschi, 2007). 
Increased mass loss, due to the increased surface opacity, can further affect the post-main sequence evolution of the 
stars, and hence the final SN yields. In contrast, it seems unlikely that rotation alone can signifi cantly alter the 



evolut ion of Pop III massive stars, without invoking anisotropic winds and strong magnetic fields ([Ekstrom et al. 



2008bD . 



Crucial for the nucleosynthetic identification of the metal-free stellar population is the uniqueness of its chemical 
signature. A clear dependence of the Pop. Ill signature on stellar mass, or ranges of masses, provides the means 
to probe the primordial IMF. Table llll shows, for each of the mass regimes given in Table U the predicted explosion 
energy (if an explosion occurs), the ejected iron mass, which of the elements are expected to be synthesized and 
ejected, and the associated characteristic nucleosynthetic signatures, if any. The box roughly indicates the range of 
primordial stellar masses as predicted by theory/simulations. It is encouraging that the chemical signatures in col. 
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5 show significant variation with stellar mass. In principle, these signatures can tell us a great deal about the first 
stars, assuming the models are realistic. 

Our working definition of Pop III demands zero metallicity, but we here sometimes refer to work done at higher 
metallicities, \i Z — models are not available, or if it can further illustrate the physics involved. In addition, we 
here do not discuss the important problem of how nucleosynthesis proceeds in Pop III binary stars, or in those with a 
significant degree of rotation, and refer the reader to the speciafized literature (Yoon et al., 2008, 2012; Meader and 
Meynet, 2012). 



1. Asymptotic giant branch stars: mj < 9 



Detailed yields of intermediate-mass stars are in general very difficult to esti mate due to several poorly cons trained 
physical processes such as overshooting, mixing and dredge- up, and mass-loss ( Campbell and Lattanziol . boOSi ) . Some 
things, however, we do know. Intermediate- mass stars cannot be a source of e.g., iron-peak elements, as they never go 
beyond the helium burning stage. On the other hand, lighter elements, in particular ^^C, ^^N, and ^^O, are produced 
efficiently by the CNO cycle and 3a reactions. These are brought up to the surface by deep mixing processes, 
specifically the TDU, and the He-FDDM (also named dual core flash, see Campbell and Lattanzio, 200 8) operating 
in lower mass stars. Generally, the C/0 raticH exceeds unity. Although metal-free stars ^ 2 M.q ( Siess et al. 



I2OO2I) can efficiently convert ^^C to "'^'^N through the CN-cycle at the base of the convective envelope, which is deep 
enough to penetrate into the H-burning shell (so called hot bottom burning), the (C -I- N)/0 remains above unity. 
As a consequence of the presence of carbon (through so called "carbon injections") in the H-burning shell, the Pop. 
Ill low- and intermediate-mass stars continue to evolve as normal AGB stars with thermal pulses similar to their 
metal-rich counterparts (Siess et al., 2002, however see e.g., Fujimoto et al., 1984; Chieffi and Tornambe, 1984). 
During these pulses, ^'^Na and ^^'^^Mg are produced, as well as neon and some aluminium. Interestingly, s-process 
nucleosynthesis may occur also in primordial stars, despite the absence of heavy seed nuclei, i.e., ^^Fe. Instead, the 
isotopes of C, N, O, F, and Ne act as seed nuclei. The relatively large neutron production in these stars result in high 
neu tron-to-seed nuclei ratios which ultimately favours the synthesis of the heaviest s-process elements, like Pb and 



Bi (jBusso et all . l2001t iGorielv and Siessl . l20Qlt ISiess and GorieM 120031 ). If so, an overproduction of Pb relative to 
Ba, for example, is expected, contrary to what occurs in more metal-rich AGB sta rs. Furthermore, stars exper iencing 
hot bottom burning can produce ''Li through the Cameron- Fowler mechanism ( Cameron and Fowleii . Il971 ). This 



mechanism becomes active when ''Be, produced by a capture on '^He in the interior, occasionally is brought up to 
cooler regions by deep convection, captures an electron and forms ''Li. Whether this Li, however, survives long enough 
to be expell ed, depends critically on the unknown amount of mass-loss at the time when the abundance of Li is high at 
the surface ( Siess et aZI. 2002). This pr oduction bra nch may partly be resp onsible for the rise from the Spite plateau 
at [Fe/H] > — 1 ( TravagHo et a/.l . [2001 ). In contrast, Sbordone et al. ( 2010( ) recently called attention to the apparent 
presence of a depression from the Spite plateau at the lowest metallicities, implying that Li have been destroyed (i.e., 
assuming that the plateau indicates the original Li abundance). However, whether this destruction is due to internal 
stellar evolution within the observed sample or due to ejection and mixing of processed material from a primordial 
stellar population (see, e.g., Piau et al., 2006) is not yet known. Additional constraints and challenges are provided 
by the ^Li abundance in metal-poor stars; we refer the reader to Asplund et al, 2006 for a comprehensive discussion. 



2. Electron capture supernovae: 9 < mj Mq < 10 



The electron capture, or ONeMg core collapse SNe, have been largely neglected in models of Galactic chemical 
evolution. Few yield calculations are found in the literature, presumably because these stars synthesize and eject only 
small amounts of metals, in particular Fe (i.e., ^^Ni). Hence, the chemical evolution for the majority of the main 
elements is, in general, unaffected by the presence of electron capture SNe. In a metal-free or extremely metal-poor 
environment, however, the impact of the electron capture SNe may be significantly higher, at least in a fraction of the 
ISM volu me. This fraction depends on the ratio of elec tron capture SNe to Fe core-collapse SNe. Various estimates of 
this ratio ( Poelarends et al . 2008; Wanaio et al. . 20091 ). based on solar metallicity models and a Salpeter IMF, quote 
a number between ~ 0.04 and 0.4. Despite this relatively low ratio, it may turn out to be of importance to account 
for electron capture SNe in low-metallicity environments, in particular for certain elements (see Sec. IVI.Cp . 



Here, the ratio C/0 = nc/no denotes the absolute number density ratio o f C an d O, respectively. Hence, a ratio of C/O = 1 is 
equivalent to [C/O] = 0.26, given the new solar abundances by I Asplund et~al. 1 1 I2OO9I) . 
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The amount of ^^Ni ejected by electron capture SNe is very small, less than 0.01 A4q. IWanaio et al. (l2009l) 
estimated a final ^^Fe yield (after decay) of 0.002 — 0.004 Mq based on a 8.8 Mq, solar metallicity model. Even 
smaller amounts of carbon and a-group elements are produced, which result s in [C, a/Fe] <j; 0. This goes along 
with small light-to-heavy a-group element ratios. Interestingly, the model bv lWanaio et al. I (l2009l) predicts a large 
production of heavy iron-peak elements, such as Ni and Zn. It is argued that the electron capture SNe could be 
the dominant source of Zn in the Galaxy. Furthermore, the elements ''"Ge and ^"Zr are found to be produced in 
abundance. Indeed, the electron capture SNe could be a potential site for the nucleosynthesis of the rapid neutron 
capture (i.e., the r-process) elements. In particular, the heaviest r-nuclei (with mass numbers A > 130) have been 
proposed to form in electron capture SNe via the prompt explosion channel (Wheeler et ai, 1998; Wanajo et at, 
2003, however, see Kitaura et al., 2006) or in the shocked surface layers of the ONeMg core (Ning et al, 2007, but 
see Janka et al., 2008). The exact site of the r-process remains unknown. It also remains unclear whether there is 
more than one r-process site (Wasserburg and Qian, 2000; Qi an, 2003; Sneden e t aL, 2008; see also Sec. IIII.B.3p . as 



has been argued for the Solar System meteoritic abundances ( Wasserburg et al . 199a) 



3. Core collapse supernovae: 10 < mj Mq < 40 



The classical core collapse SNe are found in the mass range 10 < m/M.Q < 25. They are believed to be re- 
sponsible for the bulk of the Galactic inventory of intermediate-mass elements, like oxygen. One of the most 
characteristic chemical signatures of core collapse SNe, including Pop. Ill, is the enhancement of a-elements rel- 
ative to Fe ICh icffi an d Limongi l2004t iHeger and Wooslevl . I2OIOI : iKobavashi et all . I2OO6I : lUmeda and Nomotol . l2002t 
Wooslev and Weaverl . ll995| l. Due to the small neutron excess, the underproduction of odd-Z elements, like Na and Al, 
is pronounced as compared to the solar metallicity case. Neutron-rich isotopes are also produced to a lesser extent. 
Furthermore, explosive nucleosynthesis can contribute more to the yield as a consequence of t he lower mass-loss in 
Pop. Ill stars (cf. Ekstrom et al., 2008b). If so, the production of iron-peak elements is enhanced (jlleger and Wooslevl . 
2OIOI ) as compared to more metal rich core collapse SNe. The nucleosynthe sis of nitrogen is generally n ot realized in 
non-rotating stars (but see note by Heger and Woosley, 2010). As shown bv lMeynet and Maeder (2002), rotation can 
radically boost the production of ^^N in extremely metal-poor massive stars as a result of rotation-induced mixing. 
To some extent, excess ^^C (see Fig . [5]) is also synthesized by this mechanism, as well as -'^'^C, ^^'^®0, and ^^Ne 



( Ekstrom et all l2008at iHirschil . 1200 



m 



The more massive core collapse SNe (25 < m/M.Q < 40) are increasingly affected by the deeper gravitational 
potential of the star as material falls back onto the forming black hole. The two identified branches have two 
distinctly different nucleosynthetic signatures. The weak SNe are generally characterized by very high C,0/Fe-peak 
ratios and slightly higher a/Fe-peak ratios as compared to core collapse SNe in the 10 — 25 Mq mass range. The 
deficiency of Fe-peak ele ments in the ejecta r e sults fr om fallback of the innermos t layer s. Using the fallback/mixing 
algorithm developed by IPodsiadlowski et al. (2002), Umeda and Nomotol ( 20021 2005f) proposed a mixing- fallback 
scenario in which the complete and incomplete Si-burning layers are presumed to be mixed during the explosion by 
Rayleigh- Taylor instabilities prior to the fallback, in order to obtain Zn, Co/Fe ratios in accordance with observations. 
In weak SNe, the complete Si-burning region does not reach high enough temperatures to produce Zn and Co to 
the required level of abundance. In HNe, however, higher temperatures are realised and both the complete and 
incomplete Si-burning regions are shifted outwards in mass which pror note the synthes i s of V , Go, Cu, and Zn, as 
well as '^'^Ti and '**Cr which decay into '^''Ca and "^^Ti, respectively ( Tominaga et al . 2007t lUmeda and Nomotol . 
2OO4). On the other hand, Mn, Gr/Fe ratios are low in HNe, as compared to normal SNe (Fig. [5]). Moreover, as 
the oxygen-burning region increases w ith increasin g explosion energy, oxygen-burning products like Si and S become 
more abundant relative to O and G (lUmeda et a t. 200 2^). Howev er, the observational need for core collapse SNe 
with high explosion energies is debated (Heger and Wooslev . 2010f) . Gore collapse SNe with fallback may also be a 
possible site for the r-process (iFrver et a/.l . l2006h . When material falls back onto the neutron star, it gets strongly 
shock-heated. Some fraction of this material gains enough energy to escape the potential well of the neutron star 
and as it bubbles up through the accreting gas, it cools down quickly to a few x 10^ K by adiabatic expansion. When 
the expansion slo ws down, the shock-heated ejecta enters a simmering phase, in which nucleon capture may occur. 
Frver et al. I (|2006h obtain a ratio of [Ba/Eu ]= —0.2 in the total ejecta, which is close to the r-process ratio of [Ba/Eu] 



that the solar r-process value coincides with that in the metal-poor stars (|Arlandini et al. 



-0.6 inferred from observations (jBarklem et all 120051 : iFranqois et all 120071) of extreme l y met al-poor stars. Note 

. I1999I ). 



other potential r-process sites asso ciated with core collapse SNe, such as i ieutrino-driven winds (lOian and Wooslev . 
Il996l : IWooslev and Hoffrnan . 1992) and neutron star- neutron star mergers ( Freiburghaus et al. . 1999f ). It is. however, 
beyond the scope of this review to discuss the precise conditions for n-capture nucleosynthesis and to assess which 
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Initial stellar mass (solar masses) 



FIG. 5 A selection of SN yields {Z = 0) as a function of initial stellar mass. Light colors denote light and intermediate mass 
elements (C, Na, and Mg) while dark color s denote heavy elements (Mn, Fe, a nd Zn). T he specific c olor coding is displayed 
in the figure. Yields by the Japanese group (|Nomoto et ffl/.l [20061 : lOhkubo et aZ.! . [2006: Ume da and No moto. 2002) are denoted 
by squares, yields by Heger and Woosley (2002, 20 08) are denoted by dots, yields bv .Chiefli and Limongii (j2004|) by diamonds, 
and yields (only carbon) bv lEkstrom et all (|2008al ) by triangles. Yields of core collapse SNe with masses > 40 A4q are shaded 
and connected with dash-dotted lines. Yields of HNe are connected by dashed lines. Note the sometimes large discrepancy 
between different models, which can be regarded as a measure of the uncertainties associated with these calculations. 



astrophysical site(s) ma y predominantl y be resp onsible for the in ventory of r-process elements in the Galaxy. This 
can be found elsewhere ( Arnould et al. |,|l007; Sn eden et adl2008[) . 

We emphasise that there exist very different approaches to mixing and fallback exist in the literature. The Japanese 
school parameterise both the mixing and fallback, making specific mass cuts and assuming the mixing is homogenous, 
and the fallback is complete within a given mass coordinate. Thus one can essentially develop a model that is 'made 
to order' which provides limited ins ight into the physics. T h ese m odels do not specify wh at does the mixing or fixes 
the amount of fallback. In contrast. I Joggerst et all ( 2010al lb. 20091 ) take the ID models of Heger and Wooslev ( 2010l ) 
and map these onto a higher dimensionality grid. They follow the Ray leigh- Taylor instability during the explosion 
(providing mixing) and trace the fall back of material, with the most recent models including the effects of stellar 
rotation. 



4. Pair-instability supernovae: 140 < m/AI© < 260 



In contrast to the core collapse SNe, the explosion mechanism for PISNe is well understood and there is no issue with 
fallback (Heger and Woosley, 2002; Fowler and Hoyle, 1964). As a consequence, following the remark above, PISN 



16 



I5Mg, Z = 0, Ej, = l, "Ni = 0.07, M^„,{i)= 1 .76, U^^(t) = 2.02U^ 200Mg, Z=0, "Ni=7.2SM3 




Mass Number Moss M umber 

FIG. 6 The abundance pattern in the ejecta (after radioactive decay) of a 15 A4q (left) and a 200 Mq (right) model star 
with initial zero metallicity, normalized to the solar abundance. Isotopes of even-Z (black) and odd-Z (green) elements 
are connected by solid and dashed lines, respectively. Note that the two patterns are displayed on different scale. Here, Si and 
S are clearly over-abundant with respect to C and O for the 200 A4q star, while, e.g., the Zn/Fe ratio is significantly lower. 
Adopted from lUmeda and Nomotd l|2002l '). 



yields may therefore be considered more reliable than the yields of core coll apse SNe, apart from the e ffects of rotation 
and magnetic fields. Due to the very small neutron excess in their interiors (jHeger and Wooslevll2002l ). stars exploding 
as PISNe should exhibit a pronounced odd-even effect in which the ejecta show particularly low abundance ratios 
of odd- compared to even-Z elements. Interestingly, this effect is less distinct in the models by Umeda and Nomoto 
(2002, see Fig. Furthermore, the lack of excess neutrons in ad dition to less-rapid expansion timescales during 



the e xplosion, inhibit the production of rapid n-capture elements (jHeger and Wooslevl l2002t lUmeda and Nomoto . 
12002). Likewise, the production of s-process elements is inhibited in PISNe. As a consequence of the e xtended 
oxygen-burning region, elements such as Si and S are significantly over-produced (|Umeda and Nomotd . 120021 ) leading 
to highly super-solar ratios of, e.g., [Si/0] and [S/C]. The combination of these and other abundance characteristics, 
such as low Mg/Ca ratios (true over the majority of the PISN mass range) and very low ratios of Co, Zn/Fe make up 
a chemical signature unique to the PISNe (Fig. 



5. Supermassive core collapse supernovae: 260 < raj < 10^ 



Rotating stars in the mass range 260 "S, m/ M-it, 10^ ^^7. escap e direct black hole formation if an accretion disk 
is formed around the central remnant ([Shibata and Shapiro! . 120021) . If so, freshly synthesized materi al is allowed to 
be eje cted in jets towards the polar direction, after core collapse. In their 500 — 1000 Mq model stars. lOhkubo et al. 



(|2006l ) showed that > 20% of the mass undergoes explosive Si-burning with large amounts of Fe-peak elements being 
produced, in particular elements heavier than Fe, such as Co and Zn. On the other hand, relatively little intermediate- 
mass elements are synthesized and [C, a/Fe] <C are realized in the ejecta. In a sense, the yields of 260 — 10^ M.q core 
collapse SNe anticorrelates with those of normal core collapse SNe, especially the weak SNe. Overall, [X/Fe] decreases 
with heavier elements X for normal (weak) SNe while the oppo site is true for 260 — 10^ Mq SNe. Moreover, the 
odd-even effect in the ejecta of the 500 — 1000 Mq SNe studied bv lOhkubo et al. ( 2006[ ) is pronounced as compared to 
their lighter cousins. We will further discuss the contribution from various SNe to the individual chemical abundance 
patterns of the most metal-poor Galactic halo stars in Sec. IIVI 



C. Metal transport and mixing 



The implications of ISM mixing on possible anomalies in stellar chemistry have been recognized for a long time 
( Norris and BesseUL 1978 : Searle . 1977 ). So how did the first metals mix into the IGM, and into the gas clouds that 
would eventually give rise to the first Pop II stars? This is an extremely complex problem, and we are only now 
beginning to get a handle on it. Again, progress relies on advances in theoretical modeling and numerical simulation 
techniques in combination with accurate and precise observations. The early metal dispersal problem can be separated 
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into large-scale, coarse-grain mixing, and subsequent fine-grain mixing on small scales. We discuss these two aspects 
in turn. 



1. Early transport: The supernova explosion 



It is likely that Pop III.l stars will encounter different fates, according to their mass (see Table |T|, including both 
core-collapse and pair-instability SNe. We will here illustrate the feedback from such explosions by discussing the 
evolution of a PISN remnant in greater detail. The overall hydrodynamics, however, would be very similar for the 
case of a hypernova, where metal yields are much lower, but explosion energies can be similarly large (Umcda and 
Nomoto, 2002). 

After the death of a Pop III.l star as a PISN, provided that t he progenitor had the appropriate mass (see Table H), 



the resulting blast wave goes through four evolutionary stages (jGreif et aLl . 120071 ). At very early times, the remnant 
enters the free-expansion (FE) phase and propagates nearly unhindered into the surrounding medium. It expands 
with a constant velocity Uoj given by v'^- = 2_Es„/Mej, where M^j = M*, as is appropriate for the complete disruption 



encountered in a PISN. The duration of the FE phase is given by Uc = ^fc/woj, or: 



2-Esi 



where rfc is the radius at which the swept-up mass equals the mass of the original ejecta, i.e: 



V 47rmH»T-H / 



1/3 



(1) 



(2) 



where X = 0.76 is the primordial mass fraction of hydrogen, mn the proton mass, and nu the number density of 
hydrogen nuclei. After tic, the ine rtia of the swept-up mass b ecomes important, and the shock undergoes a transition 
to the Sedov-Taylor (ST) phase (jSedovl . Il946l: iTavloii Il950l ). Typical values are rfc < 20 pc and < lO'^ yr (see 
Figure [71). 

The remnant next enters its ST phase. The photoheating prior to the explosion leads to a much reduced circ umstellar 
density, with an almost flat profile of nn — 0.5 cm~'^ in the vicinity of the progenitor star (IGreif et aZ.1 . 120071 ). The SN 
blast wave approaches after about 10^ yr, catching up with the previously established photoheating shock, where 

the outlying density profile becomes isothermal. Subsequently, the ST phase ends and the remnant undergoes a second 
transition. At the high temperatures behind the shock, cooling is due to H and He collisional ionization, excitation 
and recombination processes, bremsstrahlung and the inverse Compton (IC) effect. IC cooling and bremsstrahlung 
are important at high temperatures, whereas H/He line cooling becomes significant below 10^ K. Radiative losses 
affect the energetics of the SN remnant after about 10^ yr. At this point the shocked gas separates into a hot, interior 
bubble with temperatures above 10^ K, and a dense shell at 10"* K, bounded by a hig h pressure gradient, thus leading 



< 



10 Myr (|Greif et all |2007|). With energy conservation no 



to a multi-phase structure, which remains intact for 
longer valid, the ST phase ends. 

After the ST phase ends, the pressurized, interior bubble drives an extremely dense shell, often called a pressure- 
driven snowplow (PDS). This new phase can be described with an equation of motion: 



d{MswVsh) 
dt 



(3) 



where Pb is the pressure of the hot, interior bubble. Since Mgw oc rsh in an r' 



density profile, and Pb r^^ in the 



2/5 

adiabatically expanding interior, one can solve the above equation with a power-law of the form rgh oc i^,^ , yielding 
the same scaling as the ST solution. When does the PDS phase end? The pressure directly behind the shock after 
10^ yr can be estimated to be Pb/^B ~ 3 x 10^ K cm~^ ( Greif et al . 2007). With Pb oc r^^^ and Vgh oc tf^^, one finds 



Pb cx t^^, implying that after roughly 1 Myr the interior pressure has dropped to Ph/k^ ~ 3 x 10^ K cm^'^. At this 
point pressure equilibrium between the hot interior and the dense shell has been established, and the shock is driven 
solely by its accumulated inertia. With the pressure gradient no longer dominant, the SN remnant is driven solely 
by the accumulated inertia of the dense shell, and becomes a momentum-conserving snowplow (MCS). The shock 
position as a function of time can be obtained by solving Eq. ^ in the absence of a pressure term. Since the shock 
has not yet propagated beyond the surrounding r^^ density profile, this yields an initial scaling of rgh oc t].(^ . At later 
times, the shock finally leaves the host halo and encounters neighbouring minihalos and underdense voids. Figure [7] 
summarizes the 4-stage remnant evolution. 
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FIG. 7 Evolution of a Pop III PISN remnant, from 2: ~ 20 when the Pop III.l star explodes, to z ~ 12 when the shock finally 
stalls (adopted from Greif et at, 2007). The remnant lasts about 200 Myr. The black dots indicate the mass- weighted mean 
shock radius according to a simulation, while the dashed line shows the analytic solution. The final mass-weighted mean shock 
radius is 2.5 kpc. The shaded region shows the radial dispersion of the shock. Its increase towards late times indicates the 
growing deviations from spherical symmetry when the blast wave runs into the cosmic web, surrounding the host minihalo. 



Mixing of enriched material with gas in existing star- forming minihalos is generally inefficient, indicating that the 
dispersal of metals can only occur via expulsion into the IGM. The bulk of the SN s hock propagates i nto the voids 



surrounding the host halo, and chemical enrichment proceeds via the same channel (|Pieri et all 120071) . During the 



first 10 Myr, the initial stellar ejecta expand adiabatically and apparently do not mix with the surrounding material. 
In reality, however, the high electron mean fre e path behin d the shock leads to heat conduction and gas from the dense 
shell evaporates into the hot, interior bubble ( Gull . [iQTSi) . The onset of Ray leigh- Taylor (RT) and Kelvin-Helmholtz 
(KH) instabilities is predicted to mix the metals efficiently with primordial material evaporate d from t he dense she ll, 
reducing the metallicity of the interior by a factor of a few to at most one order of magnitude ( Madau et ai, 200ll) . 

Due to inefficient cooling, the evolution of the metal-enriched, interior bubble is governed by adiabatic expansion, 
and preferentially propagates into the cavities created by the shock. Once the shock leaves the host halo and becomes 
highly anisotropic, the interior adopts the same behavior and expands into the voids surrounding the host halo in the 
shape of an 'hourglass', with a maximum extent similar to the final mass- weighted mean shock radius. The interior 
becomes substantially mixed with the initial stellar ejecta. When the shock finally stalls, the interior bubble is in 
pressure equilibrium with its surroundings, but stays confined within the dense shell. To investigate the importance 
of RT instabilities in this configuration, one can estimate the mixing length Art for large density contrasts between 
two media according to 



Art c± 2Trgt. 



sh ' 



(4) 



where g is the gravitational acceleration of the host halo. Mixing between the dense shell and the interior bubble 
takes place on scales ^ 10 pc in the course of a few 10 Myr. Such mixing, however, is generally inefficie nt, and that 



Greif et al. 



much larger potential wells must be assembled to recollect and mix all components of the shocked gas 
2008 ) ■ Turbulence arising in the virialization of the first galaxies could be an agent for this process 
2008 ; Wise and Abel , 20081) . Recent simulations of first galaxy formation indicate that cold accretion streams deeply 



Greif et 
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FIG. 8 Turbulent flows inside the first galaxies. Specifically, the host halo here has a virial mass of ~ 10 A4q, and collapses at 
z ~ 10. Shown is the Mach number in the central few kpc of the galaxy. The dashed line denotes the virial radius at ~ 1 kpc 
from the center. The Mach number approaches unity at the virial shock, where gas accreted from the IGM is heated to the 
virial temperature. Inflows of cold gas along filaments are supersonic and generate a high level of turbulence at the center of 
the halo, where typical Mach numbers are between 1 and 5. From Greif et al. (2008). 



penetrate into the center of the corresponding potential well, giving rise to high Mach number, turbulent flows (see 
Fig.©. 

The minimum mass necessary to recollect the hot and underdense post-shock gas residing at T ~ 10'^ K and 
nn ^ lO^'^'^ cni^-^ is at least M^h- ^ 10^ ■Mq. The final topology of metal enrichment could be highly inhomogeneous, 
with pockets of highly enriched material on the one hand, and regions with a largely primordial composition on the 
other hand. Recent simulations have reached a level of realism, in terms of taking into account the radiative feedback 
from individual Pop III stars, and modeling the transport of metals from the first SNe, to give us a much better 
understanding of early enrichment, at least on small scales (1 Mpc comoving). In particular, they provide us with 
hints on the l evel and homogeneity of met als inside the first galaxies, where second generation star formation will occur 
( Greif et al . 2010l : Wise and Abell . l2008l ). It is interesting that even a single SN in one of the progenitor minihalos 
already leads to a significant metallicity inside the first galaxies (see Fig. [9]). Specifically, if such a minihalo hosted 
a PISN, the resulting average enrichment would already be ~ 10"'^ Zq (see Fig. [TU)). If this is true, the J WST is 
unlikely to see any truly metal-free galaxies, even in very deep fields. 

To further address the fine structure of metal enrichment, we now discuss the physics of fine-grain mixing. 



2. Late transport: Turbulence 

Eventually, the characteristics of the SN blast wave as a dynamical system will be lost as it is overtaken by the 
random motions of the surrounding medium. These gas motions are maintained by constant energy injections into 
the medium, e.g., by stellar winds, SN explosions, and by the conversion of potential to kinetic energy during the 
collapse of star-forming clouds and are generally turbulent in nature. At early cosmic times in particular, gas motions 
were also generated by dark matter halo interaction. 
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FIG. 9 Chemical enrichment in the wake of the first supernovae. Shown are the gas density, temperature, and metallicity in 
the central 100 kpc (comoving) of the simulation box {left to right). From top to bottom: time series showing the simulation 15, 
100, and 300 Myr after the SN explosion. The inlays show the central 10 kpc (comoving) of the emerging galaxy. The metals 
are initially distributed by the bulk motion of the SN remnant, and later by turbulent motions induced by photoheating from 
other stars and the virialization of the galaxy itself. The gas within the newly formed galaxy is already substantially enriched. 
From Greif et al. (2010). 



There is no existing general theory for hydrodynamic turbulence^, particularly if the turbulent flow is inhomogeneous 
and non-isotropic. However, from a statistical viewpoint, turbulence may be seen as a random walk process and the 
transport of matter in a turbulent medium can be described as a diffusion process (jTavloil 119211) , similar to molecular 
diffusion. The derivation outlined bel ow of the relation between time and the mean displacement of a fluid element 
in a turbulent medium follows that of Choudhuri (Il998h . 

Assume that the three dimensional displacement, x, of a fluid element after time t' = t is 



x(t) = / VL(<')d<', 



(5) 



where vl(<') is the (Lagrangian) fluid velocity at the position of the element at time t' . The mean square displacement 
is thus given by 



t t 



{x'it)) = (x(i)-x(i)) =JJ (vL(i')-VL(t"))di'di" 





(6) 



This follows since, 



See: |http:// www.claymath.org/millenniuin/Navier_Stokes_Equations/ 
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lO""* 10~^ 0.01 0.1 1 10 100 

FIG. 10 Metallicity-density relation in the early Universe. Mass- weighted, average metallicity as a function of gas density, 
shown 15 {solid line), 100 {dotted line), and 300 Myr {dashed line) after the SN explosion. At early times, a distinct correlation 
is evident: underdense region are highly enriched, while overdense region remain largely pristine. Once the potential well of 
the galaxy assembles, metal-rich gas becomes dense and the relation flattens. From Greif et al. (2010). 



t t 



( J VLdi'.y VLCli") ^{JJ VL-VLd<'di") = JJ (vL-VL)di'di". (7) 







The function {vhit') ■'Vh{t")) is the velocity correlation function of a fluid element measured at times t' and t" . 
Assuming that the turbulence does not change character with time (i.e., steady-state), this correlation function may 
be expressed as a symmetric function, TZ, of t" — t' alone. 



(vL(^')•VL(^")) = («')7^(^"-^'), (8) 

where {v^) denotes the mean square velocity. At t" — t' ,we have maximal correlation and TZ = 1, while for \t" — t' \ > 0, 
the velocities associated with the fluid element at times t' and t" become increasingly uncorrelated and TZ approaches 
zero as I i" — t' |— > c». From this definition of TZ, it is possible to define a turbulent correlation time, Tcorr, such that 

-t-oo 

7^(^')d^'. (9) 



Now, for times t ^ Tcorr, the limits of integration of the inner integral in Eq. ([6]) can be expanded to ] — 00, +00 [ 
Making use of Eq. (jH]) and the fact that TZ is symmetric, the mean square displacement may be written as 



t -t-00 t 

„2/.\\ _ I I /„,2\T-,/.// _ o_ /„.2\ / jj-" 



{x\t)) = J J {v')TZ{t" - t')dt'dt" = 2Teorr(«') j dt" = 6Aurbi, (10) 
-oo 

where 2Tcorr('^^) = SDturb- Equation (jlOp describes, in a statistical sense, the passive transport of matter in a turbulent 
medium. The coefficient I?turb is identified as the turbulent diffusion coefficient and is given by 

+ CXO 

Aurb = Tcorr (^'^ ) /3 = i(w') J TZ{t')dt' = (11) 
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where /corr = UraisTcon- denotes the turbulent correlation length and Wrms = \/ {v"^) is the root mean square velocity. 
The above expressions for the mean square diplacement and -Dtuib corresponds exactly to those of molecular diffusion. 

The continuous macroscopic "stirring" of the gaseous medium is paramount to the transport and mixing of metals 
over large distances. As compared to the microscopic, molecular diffusion coefficient, the turbulent diffusion coefficient 
is typically several orders of magnitude larger (Bateman and Larson, 1993; Oey, 2003; cf. Pan and Scalo, 2007). Hence, 
in the presence of turbulence, the mixing of gas is greatly amplified. However, it should be kept in mind that turbulence 
still operates on a macroscopic level and only molecular diffusion is able to drive the mixing on the smallest scales. 
Turbulence merely increases the area-to- volume ratio so that molecular diffusion effectively becomes more efficient. 
The situation is similar to that of the human lung. Without the hundreds of millions of alveoli which drastically 
increase the contact area to the fine network of blood-vessels, the small volume of our lungs would make it very hard 
to breathe. 



D. Transition to Population II star formation 



Chemical enrichment by the first SNe is among the most important processes in the formation of the first galaxies. 
Efficient cooling by metal lines and thermal dust emission regulate the temperature of Pop II star-forming regions in 
the first galaxies. The concept of a 'critical metallicity' has been introduceifl to cha racterise the t r ansition of the 
star-formation mode from pr edominantly high-mass. Pop HI, to low-mass Pop II stars ( Bromm et ai . 2001; Omukai 
2OOOI : Spaans and Silk . 2005 ). Cu rrently, two competing models for the Pop III - P op II transition are discussed: (i) 
atom ic fine-structure line cooling ( Bromm and Loeb j_ 2003 ■ Santoro and ShulJ. 20061) and fii) dust- induced fragmenta - 

i-i, — Hr^i„„i, „-t „;l lonnol It~» 1 — „-t „;l loni 1 . rA™,,i,„; „-t „; lonnc . c„\. — ;j — „i „i\ lonnd. Iti„ ;i,„ „„J o™,,i,„;l lonnch 



tion (|Clark et all l2008t lOopcke a/.l . I2OIII : lOmukai e^aLl . 120051 : ISchneider et aZ.l .l2006': 'Tsuribe and OmukaiV2006') . 
Within the fine-structure model, C II and O I have been suggested as the main coolants (Bromm and Loeb, 2003), such 
that low-mass star formation can occur in gas that is enriched beyond critical abundances of [C/H] 
and [0/H]^ 



-3.5 ±0.1 

-3 ± 0.2. For the atomic phase, one can robustly assume that carbon is singly ionized due to the 



crit 



strong UV background just below the Lyman limit predicted for the early stages of reionization. The dust-cooling 
model ( Omukai et ali 20051) . on th e other hand, predicts critical abundances that are typically smaller by a factor of 
10 — 100. Santoro and Shull ( 2006[ ) have investigated the additional cooling provided by Sill and Fell fine-structure 
lines, finding that their inclusion only leads to small corrections. Dust cooling would only become important at higher 
densit ies, without being able to influence the evolution of primordial gas close to the characteristic, or 'loitering'. 



state (|Bromm and Larsonl 120041) with a temperature and density: Tchar — 200 K and richar — 10 cm Ultimately, 



this theoretical debate has to be decided empirically by probing the abundance pattern of extremely metal-poor, 
second-generation stars. 



To facilitate any such comparison with stellar abundance data, it is convenient to introduce a criterion (jFrebel et 



20071 ) for low-mass star formation that combines cooling due to C II and O I. An intuitive way to arrive at the so-called 
'Frebel criterion' (following the detailed calculation in Bromm and Loeb, 2003) is to consider the balance between 
fine-structure line cooling and adiabatic compressional heating: 



At, 



Acii + Aoi >r 



ad 



where all terms have to be evaluated at nchar, ^chai - Heating due to adiabatic compression is given by 

^B^cha 



ad 



1.5n 



char " 



tff 



2 X 10 ■''^ erg s"^ cm^^ 



(12) 



(13) 



where tg ^ b x 10^ yr is the free-faU time at the characteristic state. The coohng terms are (|Stahler and PallaL l2005l) 



and 



Aoi ^ 2 X 10-20 erg s'^ cm-^ / ('2^ 



Acii 6 X 10-2" erg s-^ cm-^ ( !^ ) / ( !^ 



(14) 



(15) 



There may exist a critical metallicity which distinguishes the formation of individual stars from star clusters. This may relate to the 
critical column density that a collapsing cloud must reach through effective cooling before the cluster can form, and may explain therefore 
why all globular clusters regardless of age have a metallicity above [Fe/H] ~ —2.5. 
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FIG. 11 Testing theories for the critical metallicity with stellar archaeology. Shown is the transition discriminant, Dtrans, for 
metal-poor stars collected from the literature as a function of Fe abundance. Top panel: Galactic halo stars. Bottom panel: 
Stars in dSph galaxies and globular clusters. G indicates giants, SG subgiants. The critical limit is marked with a dashed line. 
The dotted lines refer to the uncertainty. References for the data can be found in the original journal paper. From Frebel et 
al. (2007). 



At the loitering state, C II is already in local thermodynamic equilibrium (LTE), leading to an approximate scaling 
( Stabler and Pallal 12005 ^) of tbe cooling rate: Acii oc richar exp (— 92 K/Tchar), whereas 01 is still in non-LTE, leading 



to: Aoi cx n^jjg^^ exp (— 230 K/Tchar)- Here, LTE refers to the bigh-density case, where collisions dominate tbe atomic 
level populations, wbereas non-LTE denotes the opposite situation, where radiative transitions dominate. One can 
combine these expressions as follows: 

10[C/H] + 0.3 X loi^/^] > 10-3-5 , (16) 

leading to the 'transition discriminant': 

Aran. = logio (lO^/H] + 0.3 X 10^^ /^'\ 

sucb that low-mass star formation requires -Dtrans > -Dtrans,crit — — 3.5±0.2. Using this new quantity has tbe important 
advantage that one can now accommodate an inbomogeneous data set of stellar abundances, including stars where 
both carbon and oxygen abundances are known, or even only one of them. 

However, tbis critical gas metallicity is still poorly determined, possibly depending on environmental p arameters, 
sucb as the presence of a soft UV (Lyman- Werner) radiation background ( Safranek-Shrader et all uOlut ). It is not 
even clear if there exists sucb a sharp transition. Since the enrichment from even a single PISN by a very massive 
Pop HI star likely leads to metallicities of Z > 10~^ Zq (see, e.g., Greif et al., 2007, 2010; Wise and Abel, 2008, 
see also Karlsson et al., 2008), in excess of any predicted value for the critical metallicity, these arguments may be 
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somewhat academic. The characteristic mass of prestellar gas clumps is hkely determined by a number of physical 
processes (e.g., turbulence and, possibly, dynamo-amplified primordial ma gnetic fields) other th an radiative cooling. 
The overall effect of gas metallicity on star formation may well be limited ( Jappsen et all |2009() . 

Star formation in the first galaxies may be influenced by yet another fa ctor: the CM B likely played an important 
role in setting the minimum gas temperature in the high redshift Universe (jLarson , Il998l) . The radiation temperature 
of the CMB is Tcmb = 2.728 K(l + z) such that Tcmb is greater than about 30 K at z > 10. Hence, early molecular gas 
clouds cannot cool to the characteristic temperatures of present-day star-forming regions 10 K) even if molecules 
such as CO and H2O act as efficient coolants. If the minimum gas temperature is a critical physical quantity that 
determines the characteristic masses of stars, then the transiti on from high-mass Pop III to the lower-mass Pop II 
mode may be gradual and regulated by the CMB temperature ( Schneider and Omukai . 2010l) . 

Anyhow, It is intriguing to see how well the fine-structure model, based on C and O, does in terms of accounting 
for all existing metal - poor stellar data (se e Fig. [TT|) . While the criterion has been criticised on theoretical grounds 
( Jappsen et al . 20071 : lOmukai et oil . l2005l) . the simple model matches the observed lower threshold very well, also 
in the regime below [Fe/H] ~ —4. Only one probable exception is currently known. With an Fe abundance of 
[Fe/ H] = -4.7 (ID LTE analysis), the star SDSS J102915 -I- 172927 is among the most iron-poor stars observed to 
date (jCaffau et al. Interestingly, it a ppears not to be pa rticularly enhanced in carbon, with a measured upper 

limit on the C abundance of [C/H] < —3.8 ( Caffau et all [2OII]) . Together with a presumed 0-to-Fe enhancement of 
[0/Fe] = -1-0.6 (typical for "C-normal", metal-poor stars), SDSS J102915 172927 falls into th e 'forbidden zone' in 
Fig, rm This implies the presence of a dust-cooling channel in Z ^ 10^^ Zq star- forming gas (jKlessen et all 12012 : 
ISchneider et"all . l2012l) . Since the majority of stars below [Fe/H] = —4 are enhanced in the CNO elements and fulfill 
the Dtrans Criterion, both cooling channels may in fact be viable. If so, the fraction of stars below the line-cooling 
threshold would then be a measure (although complex) of the relative importance of dust for the transition to Pop II 
star formation. Future observations should aim to constrain this fraction. Recently, evidence for carbon enhancement 
toward low [Fe/H] has also been discovered in an extremely metal-poor damped Lyman-a (DLA) system at z ~ 2.3 



(|Cooke et a/.l . 120111 ). opening up a new window into the conditions for early star formation. With a view to upcoming 
JPFST observations, it is important to formulate similar empirical tests for alternative models for the Pop HI - Pop II 
transition that exploit dust cooling, say. 



IV. TRACERS OF PRE-GALACTIC METAL ENRICHMENT - CONCEPTS 

This review seeks to identify the chemical signatures of Pop HI stars in the surface abundances of the oldest stellar 
populations. As already noted, we make a clear distinction between a chemical fingerprint and a chemical signature 
to readily distinguish between microscopic and (a series of) macroscopic processes respectively. We now discuss these 
in turn. 



A. Chemical fingerprints (microscopic processes) 



The many elements that can be measured in metal-poor stars provide evidence of physical processes that we are 
o nly now beginning to r ead ( e.g., Sned e n et al., 2008). This approach has a long tradition. The celebrated papers 
of Burbidge et al. (1957) and ICameron ( 1957 ) demonstrated how to infer the processes involved in stellar chemistry 
from detailed abundance measurements in meteorites and stars. This tradition can be traced back e ven further in the 
early interpretation of the Solar System abundances (see Suess and Urey, 1956). In a seminal paper, ICameronI (Il973l) 
showed how the meteoritic record can be separated into distinct neutron and proton capture processes. Remarkably, 
as we mention below, the fingerprints of these processes are observed directly in a substantial fraction of the most 
metal poor stars. As a general rule, we have a better understanding of fingerprints when compared to chemical 
signatures. For example, while there is still no in situ evidence of r-process nucleosynthesis in supernovae, the main 
reaction networks that give rise to different heavy elements for a given neutron flux and number density are fairly 
well understood. How the s-process ele n ients are synthesized in quiescently evolving stars is also well understood 



(|Busso et a/.l . ll999[) . After ISneden et all I (l2008l) . we emphasize that the stellar record is not complete at the present 



time. There are no elements measured from Ag to Ba in any n-capture rich stars, which will be needed to rule out 
alternative n-capture models in the most metal poor stars (see Arnould et al, 2007). Thus, chemical flngerprints are 
concerned with distinct hallmarks of a microscopic process that is broadly understood in terms of a limited set of 
physical parameters (e.g. flux, density). 
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B. Chemical signatures (macroscopic processes) 

Most of this review is concerned with identifying and interpreting chemical signatures of the first stars and sub- 
sequent generations. This is inherently more difficult. In discussing these signatures, we are largely dependent on 
the yield calculations from complex hydrodynamical simulations (Sec. IIII.Bp . The calculated yields are dependent on 
numerous factors, including stellar mass, metallicity, rotation, black-hole infall, and so forth. The numerical models 
make very detailed predictions up to atomic masses of Z « 30, which falls short of the neutron-capture elements (cf. 
Pignatari et ai, 2008), so the predictions are limited at t he present time. But these models predict distinct chemical 
signatures. For example, in intrinsically faint supernovae, Umeda and Nomotol (2003) find strong CNO enhancements 



compared to Fe due to mixing and fallback onto the black hole. In contrast, hypernova models predict strong Zn and 
Co as compared to Cr and Mn (e.g., Tominaga et ai, 2007). 

Apart from the underlying signatures of the preceding stellar generation(s), there are likely to be many macroscopic 
processes that influence the observed abundance pattern of metal-poor stars. Stars that are born into binary pairs can 
be enriched by mass transfer from the primary donor, either through slow winds or Roche-lobe overflow. A significant 
fraction of the most metal-poor Galactic halo stars show evidence of mass transfer from a companion that is no longer 
directly visible, as we discuss below (e.g., Fujimoto et ai, 2000; Suda et at, 2004). In principle, this can complicate 
our interpretation of element groups within the spectrum of a metal-poor star as the original chemical signature is 
diluted with a secondary one. Another complicating factor is mixing within the interstellar gas. The efficiency of the 
macroscopic (turbulent) mixing of the ISM can transform the underlying signature and affect the overall trend as a 
function of metallicity (i.e., [Fe/H]). 

Some signatures only reveal themselves by comparing the abundances of an ensemble of stars. One example is 
clustered data points in the abundance plane as a result of clustered star formation (which produces chemically 
uniform star clusters). As discussed in Sec. I Villi this is an extremely important phenomenon because it allows us to 
reconstruct star clusters that have dissolved in the distant past. Of no less importance, it allows us to identify stars 
that are not subject to chemical anomalies due to mixing or binary mass transfer, thereby providing us with a cleaner 
signal of the preceding generations of stars. 



C. First versus second generation stars 



In this review, we have identified the Pop III stars (see Sec. U for definition and Sec. |TT] for discussion on the 
formation of Pop III.l and Pop III. 2 stars) as being metal-free, which in the context of early star formation implies 
a metallicity Z < Zcrit (see Sec. IIILPp . If nothing else is stated, we will continue to use the terms 'primordial stars' 
and 'first stars', as synonyms for the Pop III. Chronologically, some of these designations, in particular 'first stars', 
may be slightly ambiguous since the formation of Pop III. 2 stars likely occurred as a result of feedback, e.g., from 
Pop III.l stars (Sec. HH. The Pop III.2 stars were therefore not, in a literal sense, the very first stars in the Universe, 
although they were the first stars to form locally. From a chemical point of view, however, both Pop III.l and Pop 
III. 2 stars can readily be defined as 'first stars'. 

Second generation stars are defined as stars formed out of gas only enriched by the primordial stellar generation. 
As such, they may also b e ca lled the fi rst generation Pop II stars. As primordia l star formation occurred in different 
places at different times (jBrook all 12007 : Maio et all . 2010t Tornatore et al . 2007), second generation stars may 



not necessarily belong to the oldest stellar population in the present Universe, although a fraction of the oldest stars 
should belong to the second generation. They will, however, generally be very old and probe a vast variety of star 
formation sites and initial conditions, from the Galactic bulge to dwarf galaxies. 

The identification of second generation stars is a difficult task as there is no way of labeling the heavy element 
atoms as 'synthesized in a metal-free star' or 'synthesized in a metal-enriched star' at the present time. The hope is, 
however, that the more deficient in metals a star is, the more likely it is that it belongs to a generation immediately 
following that of the first. While this should be true in a general sense, despite a significant spread in the predicted 
age-metallicity relation (e.g., Argast et ai, 2000; Tumhnson, 2010), it may not be true in particular cases, as we shall 
see (see, e.g., Sec. IVI.DI and IVIII.B|) . In any event, the identification of extremely metal-poor (EMP) stars (defined 
as stars with an iron abundance of [Fe/H] < —3, Beers and Christlieb, 2005) in the Milky Way and in other galaxies, 
is essential for probing the Pop III era. 
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FIG. 12 Simulation ol the chemistry of a Milky Way-like environment. The maps span 300 kpc on a side. Note the enhanced 
levels of [a/Fe] in the Galactic halo while the rati o is nearly s olar i n the satellite galaxies, indicating a significant contribution 
of iron from thermonuclear (Type la) SNe. From lFbnt" et al\ l|2006l) . 



V. IN SEARCH OF THE EARLIEST CHEMICAL SIGNATURES 



A. Sites of early star formation 



Essentially all galaxies, including all known dwarf galaxies, appear to contain stars older than 10 Gyr. But only a 
fraction of these are likely t o contain pristine signatur es of the first and s econd stellar generations. So where do the most 
ancient stars reside today? IScannapieco et al. ( 2006 ) and Brook et al. ( 2007 ) showed that metal- free stars, i.e., if they 
exist, should be spread over the entire Galaxy. This implies that slig htly enriched, second generation stars, sho uld be 



found more or less everywhere, and predominantly in the outer halo ([Brook et all l2007t ISalvadori et al\ . \201(X). This 



is als o where the most metal-p oor stars have been fo und in past and present surveys. Together, the HK (jBeers et 



19921 ) and the Hamburg/ESO ijWisotzki et al 1 12OOOI) surveys of metal-poor stars have discovered roughly 250 stars 
below [Fe/H] = —2 and several hundred stars below [Fe/H] = —3 in the Galactic halo ( Beers and Christliebl . 12005) . 
The Sloan Digital Sky Survey (SDSS) and Sloan Extension for Galactic Understanding and Exploration (SEGUE) have 
now raised the numbers of known stars with [Fe/H] < —2 to over 25, 000 and the number of stars with [Fe/H] < —3 to 
roughly 1, 000 (Beers, private communication). Stars with different metallicities are distributed differently in the halo. 
The outer halo of the Milky Way appears to be comprised of a stellar population significantly more metal-poor than 
its inner halo (Carollo et ai, 2007, 2010; but see Schonrich et al, 2011). According to these results, the metallicity 
distribution function (MDF) of the outer halo peaks around [Fe/H] ~ —2.2, while the MDF of t he inner halo peaks at 
[Fe/H] c± —1.6. The outer an d inner halo may also differ in terms of oblateness and net rotation (jKinman et all 12007 



MaiewskiL IToqI IZinnl . [l993l) . This suggests that the outer component was formed in a distinctly different way from 
the inner halo, probably by dissipationless merging of small dwarf galaxies at a relatively late state of the assembly 
of the Milky Way (e.g., Carollo et al, 2007, see also De Lucia and Helmi, 2008; Zolotov et al., 2009; Salvadori et 
al., 2010). As a result of the inside-out construction of the Milky Way galaxy the oldest, inost n i etal-poor stars ar e 
expected to r eside in the innermost part of the Galactic halo ([Bland-Hawthorn and Peebles! . [2006 ; lTumlinsonl . l2010[ ). 
or the bulge ([White and Springe! . I2OOOI) . However, recent developments in galactic dynamics may force us to revise 
this expected picture, as discussed below (see Sec. IV. Bp . 
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Extremely metal-poor stars are now also being found in dwarf galaxies in the Local Group. An increasing number 
of high-resolut ion spectrome t ric studies of the brightes t stars in the faintest galaxies consiste ntly find stars with 
[Fe/H] < -3 (jFeltzing et all 120091: iFrebel et all l2010bl: iNorris et all I2OIOI: ISimon et all l2010t) . Interestingly, the 



Helmi et al.. 


20061). 


Shetrone et 


al. 2001: 



In combination with observations of lo w [a/Fe] ratios in the more metal-rich member stars 



Shetrone et all I2OOII : iTolstov et al. 1. 120091 : IVenn et al. I. I2OO4I) . this casted doubt on the present-day dSphs as being 



the formation sites of Galactic halo stars. Nor did they seem to be the surviving counterparts to the now dispersed 
systems in which the Galactic halo stars once w ere formed. Howev e r, new and r efined measurements reveal th e 
existence of EMP stars in these sys tems as well (jFrebel et al ]. l2010at iKirbv et all I2OO9I : IStarkenburg et all 120101 ). 
The simulations by Font et al. ( 2006f ) depicted in Fig. [121 nicely illustrate the general understanding of the distribution 
of metals in the stellar component, in the Milky Way environment. 

At present, it is not clear in which environments the bulk of the EMP stars were formed originally. Since they are 
believed to be second generation stars, born out of gas enriched only by the first SNe, the majority of these stars must 
have formed relatively early on, possibly in so-called 'atomic cooling halos', which are small dark matter halos on the 
order of 10^ Mq (see, e.g., Greif et al, 2008). These halos were probably among the first objects able to form Pop 11 

stars (e.g., Karlsson et al, 2008; Greif et al, 2010). In fact, the majority of the stars in the newly discovered ultra- 

' — 1 1 — — ' " — — ^ ' ^ — 1 1 — — 1 1 irr 

faint dSphs ( Belokurov et aLI. 120101. 1 2007: Willman et al, 200 51: IZucker et al '. 2006a, bf), suggested to be relic galaxies 

forme d prior to reionizat ion (jBovill and Ricotti . ,2009; .Salvadori and Ferrara . ,2009,1 are very metal-poor with [Fe/H] 



< —2 (jKirbv et a/.l . 120081 ). If these galaxies should turn out to be surviving members of a population of 'first galaxies' 
connected to the atomi c cooling halos, it would g ive us a unique possibility to directly study the birth sites of second 



generation EMP stars (jFrebel and Bromml 1201 



B. Secular evolution, blurring and churning 



In the classical picture (e.g. Wielen, 1977), star clusters grad ually dissolve and diffus e into the background galaxy 
staying within a kiloparsec or so of the birth radius. But since ISellwood and Carlberd (Il984l) . we have slowly come 
to realize that spiral p ertur bations are driven by accretion events and are likely to be transient. In a seminal paper, 
Sellwood and BinnevI ([2002) showed that in the presence of such perturbations, a star can migrate over large radial 
distances. During interaction with a single steady spiral event of pattern speed Jlp, a star's energy {E) and angular 
momentum (J) change while it conserves its Jacobi integral. In the {E, J) plane, stars move along lines of constant 
Ij = E — rip J. For transient spiral perturbations, the star undertakes a random walk in the {E, J) plane, deflected 
by a series of presumably uncorrelated spiral arm events which occur on time intervals of 500 Myr or so (the impulse 
acting on individual stars lasts for a fraction of this interval) . For example, the Solar family would have experienced 
about ^ 10 of these events in the vicinity of co-rotation if only a fraction of the stars are moved in each event 
"2010bh . 



(Bland-Hawthorn et al 



Substantial variations in the angular momentum of a star are possible over its lifetime. A single spiral wave near 
co-rotation can perturb the angular momentum of a star over a broad distribution in with a tail up to 50%. The 
star is simply moved from one circular orbit to another, inwards or outwards, by up to 2 kpc or so. Conceivably, the 
same holds true for a star cluster except that most are dissolved in much less than an orbit period. 

Support for the Sellwood-Binney mechanism (churning) comes from recent N-body simulations that allow for a 
steady accretion of gas from the outer halo (Roskar et al, 2008). These relatively "cold" simulations reveal some of 
the secular processes that may be acting over the lifetime of the disk. Averaged over the entire disk, these authors 
find that the rms radial excursion at t = 10 Gyr is ((Ai?)^)'^-^ w 2 — 3 kpc. For the outermost disk stars, the 
simulations show that (Ai?) « 3 — 4 kpc, several times larger than their epicyclic radial excursion. There are related 
mechanisms that may behave even more efficiently and can even extra c t stars from the inner bulge and transport 
them across the disk ( Minchev and Famaevl . I2OIOI : iMinchev et al , 2012 ; Quillen et al , 2009t Schonrich and Binnev , 
20091) . The widespread population of supersolar metallicit y stars throughout the disk lends further support to the idea 



that circumnuclear stars can travel far out into the disk (jTrevisan et aLl . l201lh . This a ssumes that most supersolar 



stars are born in the inner regions which finds support from high redshift observations (jHamann and Ferlandl . ITooi 

2012h . Thus, it is no longer clear that the most ancient stars are highly concentrated towards the 



Savaglio et al 



centres of galaxies. 
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FIG. 13 The [a/Fe] abundance ratios in very metal-poor Galactic halo stars. The bulk of the halo sample show a clear 
enhan cement in the a-group-to-Fe ratios, in particular for Mg, Ga, Si, and Ti, as compared to the Sun. From ICavrel et al\ 



VI. CHEMICAL SIGNATURES FROM THE GALACTIC HALO 
A. General behaviour 



Detailed chemical abundances based on high-resolution {R > 30 000) spectroscopy for many ele ments (e.g. C , 
N, O, Na, Al, Mg, Ca, Si Sc, Ti, Cr, Mn, Fe Co, Ni, Zn , Sr, Y, Zr, Ba, Eu) are now availab l e dBarklem et 



2005; Bo nifacio eroH. 120091: ICavrel et ai\ . 12004 ICohen et all l2004l : [Francois et dl . 120071: iLai et all 120081: ISpite et al 
2005; .Zhang et al. . 20 111 ) for a large number of Galactic halo stars below [Fe/H] = —2. The general conclusion is 



clear: the large majority of the EMP stars were formed out of gas enriched by core collapse SNe in the mass range 
10 ^ m/MQ ^ 40 (see Sec. IIII.B[) This is concluded primarily from their enhanced [a/Fe] ratios (see Fig. [T51) . 
There are, however, a number of exceptions to this observation. The behaviour of Zin c, or rather [Zn/Fe], at the 
low est metallicities cann ot easily be explained by normal core collapse SNe (Fig. [HI). lUmeda and Nomotol (|2002| ) 
and Umeda et ^ (12002 ') showed that the production of Zn is favoured in their more energetic, hypernova models as 
a result of the extended S i-burning region (see also Umeda and Nomoto, 2005; Tominaga et al., 2007). In contrast, 
Heger and Wooslev ( 2010l ) argue against HNe as the primary nucleosynthesis site for Zn and notice that this element 
has an uncertain contribution from the s-process, inter alia. Furthermore, since only a fraction of the massive stars is 
believed to explode as HNe, a large spread in [Zn/Fe] would be expected at the lowest metallicities, i.e., if the ISM was 
poorl y mixed. Such a. scatte r is not observed, although a number of relatively low upper limits in [Zn/Fe], re ported 
on bv lBonifacio et al. ( 20091) . may possibly be suggestive of a larger scatter. In addition, Honda et al. (2011) found 
a strongly Zn-enhanced star at [Fe/H] = —3.2, that could be explained by a HN contribution. 

Carbon and nitrogen are two other elements which behave rather curiously for 'normal' stars in the metal-poor 
regime. The observed star-to-star scatter is significant for both elements (even if the population of C-enhanced stars 
are excluded), in particular for N (e.g. Spite et al., 2005). The evolution of [C/0] a t low metallicities (Fig. [14]) tends 



2004 iFabbian et ad 120091) . This has been interpreted as a 



to increase with decreasing metallic itv (lAkerman et al 

unique signature of metal- free SNe ( Akerman et all l2004li ^. However, this signature may alternatively originate from 
massive Pop II stars in which significant a mounts of C (an d N) are dredged up and ejected in an ex t ensive stellar 



massive rop it stars m wmcn signmcant ar nounts ot (an a in ) are dreagea up and ei 
wind induced by high rotational velocities (|Chiappini et al . . l2006t iFabbian et all liooa 



Mevnet et all . l2006[) . If so. 



the high C/0 ratios in the metal-poor stars are not uniquely a signature of Pop HI. 

In general, odd-Z elements li ke K, Sc, V, and Co, but also Ti, appear to be difficult to synthesize in standard 



nucleosynthesis models of SNe ( Kobavashi et al. . 20061). Irnprovem ents can be achieved by including additional nu- 



2008bD, the pi^-process (Frohhch et al. 



cleosynthesis processes, like the i^-process (jYoshida et al. 

important for the production of ^^Zn and above), or modifications of parameters such as the electron fraction 



2006, could be 
How- 
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FIG. 14 Evolution of Zink and carbon in the metal-poor, Galactic halo. The data bv lCavrel eiall (|2004H show an increase of 
the [Zn/Fe] ratio with decreasing [Fe/H] (left panel). This increase may be understood in terms of a contribution of Zn from 
HNe. The red, dashed line denotes the average slope of the data. Similarly, there is an upturn in the [C/O] ratio (right panel) 
at low metalli cities (her e measured by [0/H]). Pop III stars, or fast rotating Pop II stars, could be responsible fo r this upturn. 
Data by iFabbian et al\ ([20091') ar e denoted by black diamonds while the red, open triangles are data of DLAs bv lPettini et al\ 
(120081 ). See IFabbian et al\ l|2009l) for further references. 



ever, this discrepancy may not only be an effect of incomplete modeling of the nucleosynthesis (jHeger and Wooslev , 
2Q10I) . but could also result from systematic errors in the observations, such as non-LTE effects and /or the l ack of 



a proper 3D abundance analvsis (lAndrieysky et ali. l201 i'AsDhmdl, l2005t iBonifacio etldl . l2009t ICollet et all . l2006t 
Garcia Perez et all l200d : iLai et all 120081 ). IBonifacio et al. (2009;) identify a number of discrepancies between the 



inferred chemical abundances in EMP giants and turn-off stars. They conclude that these discrepancies are predom- 
inantly, possibly except for C, due to shortcomings in the abundance analysis. Indeed, Andrievsky et al. (2010, 
see also Takeda et al., 2009) find that a careful non-LTE analysis of Na, Al, Mg, and K tightens the abundance 
relatio ns and some agreem ent with models of chemical evolution of these elements is attained. The same goes for 



silicon tehang et al. . l201lh . It is crucial that these systematic errors are identified and reduced 



Cu riously, the signature of PISNe is not detected in the Galactic halo EMP stars (jBallero et al ]. l2006t l ICavrel et all 



20041 : iTumlinson et all 120041 ). This observation sug gests that very massive stars were ext remely rare in the early 



Universe, in contrast to some theoretical arguments (^Abel et al., '2001 Bromm et al. 



other predictions ( Karlsson et 
further in Sec. WW\ 



120091 ). but in agreement with 
al.. . ,2008; .Salvadori et aL . ,2007; iTornatore et al., . i2007[ ). This issue will be discussed 



B. Scatters, trends and mixing 



Recent observations by a number of authors (e.g., Cayrel et at, 2004; Arnonc et aL, 2005; Barklcm et at, 2005; 
Andrievsky et al., 2010; however, see Zhang et al., 2010) have revealed surprisingly small star-to-star scatters in 
the a-to-iron element ratios for very metal-poor and EMP stars, scatters which are consistent with observational 
uncertainties alone (Fig. [T5)) . This has been interpreted as evidence of a chemically well-mixed star- forming medium, 
even at very early epochs, in which the mixing timescale was extremely short. However, the star-to-star scatter in a 
given abundance ratio is only an indirect measure of the state of the early ISM/IGM. Indeed, a long star formation 
timescale produces about the same scatter as a short mixing timescale. lArnone et al. (2005) speculated that the 
longer cooling times in metal-poor gas may prevent subsequent star formation until the ejecta from whole generations 
of SNe were efficiently mixed. As a result, a small star-to-star scatter can be achieved. Very long cooling timescales 
would, however, wipe out any observed trend with metallicity, such as a decreasing [Zn/Fe] with increasing [Fe/H] 



observed for EMP giants (jCavrel et all 120041) . unless the SN yields, integrated over the IMF, show a significant 



metallicity dependence (see Fig. [15]). At least in the case of [Zn/Fe], theoretical yield calculations do not seem to 
predict such a dependence (e.g., Chieffi and Limongi, 2004; Nomoto et al., 2006). An interesting example where 
metallicity dependent yields may be able to explain observations is the decreasing [C/O] trend with increasing [0/H] 
studied by Akerman et al. , (2004; but see the comment in Sec. IVI.Ap . In their analysis, they adopted the yields by 
IChieffi and Limongil (|2004[ ). 

A viable explanation to some observed trends is a change of the IMF itself, e.g., from a top-heavy IMF to an 
IMF which more resembles the present-day one (Fig. I15|) . This route should be explored further. Interestingly, 
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FIG. 15 Predicted trends in the [Fe/Ca] - [Ca/H] plane, assuming different mixing scena rios. The light b l ue sh aded area within the blue 
solid lines denotes the 1-a dispersion in the inhomogeneous chemical evolution model bv lKarlsson et al\ l|2008l l. where SNe enrich their 
surroundings locally (although mixing will occur on a longer time scale). The negative slope is a result of the mass-dependence of the 
Ca yield where regions enriched by the least massive SNe have low [Ca/H] and high [Fe/Ca] while regions enriched by more massive SNe 
have high [Ca/H] and low [Fe/Ca], The red thick dashed line denotes the trend predicted from a homogeneous chemical evolution model 
assuming instantaneous mixing. The positive slope originates, again, from the mass dependence of the Ca yield. In this case, however, 
the most massive stars are the first to explode and enrich the (entire) system, producing a low [Fe/Ca], while the least massive SNe enrich 
the system slightly later. The thick green dash-dotted line denotes a case where the cooling time scale of the gas is long as compared 
to the mixing time scale and the ejecta of SNe of all masses have time to mix before the next generation of stars is formed. Here, the 
positive slope is a result of a metallicity dependence of the SN yields such that the IMF averaged [Fe/Ca] ratio is slightly higher at higher 
metallicities. The thin black dashed line denotes a similar case where, instead, the IMF is allowed to change with metallic ity, from a 
top-heavy to a Salpeter one (no metall icity dependent yiel ds are assumed). In all cases, the SN yields bv lNomoto et al. have been 

adopted. Data of Galactic halo giants ijCavrel et g^L [20041 ) are shown as black circles, for comparison. 



there are claims of IMF variations also in the local Universe which appear to arise in low surface brightness or low 
star-formation rate galaxies. There is now a very well established deficiency of Ha luminosity relativ e to other star 



formation indicators (e.g. UV luminosity) that is not easily explained away without modifying the IMF (jMeurer et al 
[2009) . One possibility is that ionizing photons are leaking out of low surface brightness galaxies, such that the Ha 
luminosity does not reflect the true ionizing luminosity. But the effect could herald a non-uniform IMF at the present 
or any epoch as a consequence of galaxy evolution. Possible observational effects of an evolving early IMF are discussed 
further in Sec. IVI.DI 

The presence of trends in the EMP Galactic halo may alternatively suggest that the ISM at this point was, in fact, 
fairly unmixed. In a poorly mixed, extremely metal-poor ISM, two low-mass stars can form out of gas enriched by 
two SNc of different masses producing different amounts of h eavy elements. In such a scena r io, trends may simpl y 



result from the different mass-dependences of the SN yields (jKarlsson and Gustafssonl . 120051 : iNomoto et al. , l2006h . 



However, in contrast to what is a consequence in homogeneous chemical evolution models, where the most metal-poor 
stars are enriched by the stars with the shortest lifetimes, i.e., the highest-mass stars, the most metal-poor stars will 
instead, in most cases, be enriched by the least massive SNe as they eject the least amount of metals. 

A number of elements show clear evidence of inhomogeneous chemical enri chment. Large and real scatters are found 
in the n-capture elements (e.g., Frangois et ai, 2007) as well as in C and N (jSpite et all , liooi). Although the scatter 
in the n-capture elements may result from the possibility t hat these ele i nents are only prodii ced in a small SN mass 
range as compared to light and intermediate-mass elements dBurris et aL '.'2000': Cescutti, 2008"), this would not explain 
the observed scatter in those other elements. As regards nitrogen, IChiappini et aL (,2006 ) speculated that the large 
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observed scatter may originate from variations in the initial r otational velo city of extremely metal-poor massive stars, 
which are shown to produce large variations in the N-yield (lHirschi[2007h . In general, a large star-to-star scatter is 
realised in metal-poor an d unmixed gas for elements wi t h stellar yields that depend strongly on some paramet er, like 
initial mass or rotation ( Cescutti and Chiappini . 20101 : iKarlsson and Gustafsson . 20051 : Karlsson et ai . 2008() . The 
scatter dichotomy, where some elements show evidence of an unmixed ISM while others do not, may also need to 
be understood in the context of star cluster for mation. If the majority of the E MP stars originally were formed in 



clusters, as are stars in the present-day Universe (jBland-Hawthorn et g/J . l2010a^ . extra mixing will initially occur on 



the scale of the H2 clouds out of which indivi dual clusters were for med. Although this could provide enough mixing to 



explain the very small scatter e.g., in [a/Fe] (jArnone et a/.l . l2005r ). the large scatter observed for the neutron-capture 
elements must still be accounted for. This potential 'cluster-mixing issue' requires further attention. In Sec. IVIIH 
we will come back to the importance of star clusters and the exciting prospects to use their chemical signatures as 
probes of the early star and galaxy formation. 

The issue of trends, as illustrated in Fig. (TSl and of the scatter dichotomy needs to be addressed in more detail. 
It is important to keep in mind that incomplete and/or biased samples of stars may not capture the true cosmic 
star-to-star scatter. Moreover, the effects of inhomogeneous chemical enrichment are only expected to be observed 

below [Fe/H] 3 (iKarlssonL lioiil ) . a metallicity regime which still is quite poorly sampled. Future stellar surveys 

need to go deeper, with a focus on minimizing the observational uncertainties. The signatures of chemical evolution 
in physically distinct stellar systems, such as the bulge, halo, and dwarf galaxies, can then be inter-compared and be 
used to disentangle progenitor mass-, age-, or metallicity dependences on abundance trends (Fig. I15p . In order to 
fully unveil the chemical signatures of Pop III, it is of utmost importance to gain a deeper understanding on what 
drives the interstellar gas mixing, how efficient it is in various environments (e.g., Pan and Scalo, 2007; Greif et al., 
2010), and how it couples to star formation. 



C. Outlier objects and carbon-enhanced metal-poor stars 



The Hambur g/ESO survey of met al-poor stars has discovered the th ree most iron-poor st ars known to date : 



HE 0107-5240 dGhristheb et Qi.l . l2002h . HE 1327-2326 ^^^^dl,^QQi) , and HE 0557-4840 (" Norris et oi.l . l2007h . 



All three stars are ultra-metal-poor, lying well below [Fe/H] = —4, while two of them, H E 0107 and HE 1327, are 
hyper-metal-poor ([Fe/H] < —5) according to the definition bv lBeers and Christliebl ( 20051 ). Curiously, these stars, in 
particular the two most iron-poor ones, also exhibit several abundance anomalies as compared to the bulk of Galactic 
halo EMP stars, including strong enhancements of CNO elements relative to Fe (only an upper limit of [N/Fe] < 1.0 is 
reported for HE 0557, see Norris et al, 2012), Na (HE 0107 and HE 1327), Mg (HE 1327), Al (H E 1327), and Ca (HE 
1327). Furthermore, despite HE1327-2326 being almost unevolved, it is strongly depleted in Li ( Frebel et aLl . 12008). 
Different scenarios have been put forward to explain the anomalous abundance patterns, in pa rticular the high CNO 
abundances . These scenar ios include pre-enrichment of the pri mordial cloud by a f aint S N (jlwamoto et al. I. I2OO5I: 
Umeda and Nomoto . l2003h . enrichment by mass ve stellar winds ( Mevnet oZI . l2006l 120101 ) , mass transfer from a bi- 
nary companion (e.g., Suda et al, 2004; Tumlinson, 2007a), and atmospheric dust-gas separation (jVenn and Lambertl . 
20081 ). So far, none of the stars show radial velocity variations indicative of binarity and none of them seem to be 
particularly enhanced in s-process elements, which otherwise often is the case for stars affected by mass transfer from 
an AGB co mpanion. 

Recently, ISpite et all (l201lh argu ed that the abundance pat terns of HE 0107 and HE 1327 presumably did not suffer 
from gas depletion, as suggested by I Venn and Lambert ( 2008 ). since the somewhat more metal-rich star analogue, CS 
22949-037, appears to show no enhancements in S and Zn. Enhancements in these elements would have been expected 
if the abundance pattern of CS 22949-037 had been locked to the dust condensation temperature, as it i s in stars 
affecte d by dust depletion, like the A Bootis stars. As an alternative explanation for the high C abundance, [Karlsson 
( 2006f ) explored a scenario in which the C (and N, O) originated from massive rotating stars (see also Meynet et al, 
2010) with m > 40 A^0, while the Fe originated from somewhat less massive core collapse SNe in combination with a 
period of reduced s t ar for mation after the first generation of stars had formed (see below). For a similar scenario, see 
also Timongi et all (j2003[ ). Despite the unusually high CNO and the enhancements of a few of the intermediate-mass 
elements, normal core collapse SNe in the range 15 ^ m/M.^ ^ 30 seem to be fairly successful in describing the 
abundance ratios of the heavier elements. 

It is not yet clear whether these stars tell us something fundamental about the first stars and the conditions of star 
formation in the early Universe, whether they are "outlier objects" formed under very special circumstances but not 
necessarily unique for Pop III star formation, or whether their chemical abundance patterns, at least some part of it, 
have been altered by a process subsequent to their formation, in which case the original signature would be hidden 
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FIG. 16 The predicted distribution of stars in the [C/Fe] — [Fe/H] plane. In both panels, the thick, dashed line indicates a carbon 
abundance of [C/H] = —3.5. The shaded area below this limit should, according to the 'Frebel criterion' (see Sec. IIII.DI l. contain very 
few low-mass stars and are therefore disregarded here. Symbols denote various observations of stars in the Milky way's halo (for details, 
see Karlsson, 2006). The peak of the predicted group of mega metal-poor stars (i.e., [Fe/H] < —6) is encircled. In the lower panel, the 
fiducial c arbon yield is increased by a factor of 4 for stars in the mass range 30 < ml Mq < 60, in ord er to estimate th e effect of stellar 
rotation llMevnet and MaedeJ. 120021') . Observed abundance ratios are not corrected for 3D effects. From iKarlssonI ||2006|) . 



from us. It is important, however, not to disregard them in the context of chemical evolution purely based on the 
argument that they are "chemically weird" and lie too far away from the mean. Such a decision should be made based 
on more physical grounds and detailed observations (see Sec. IVIII|) . 

Carbon plays a special role in early star formation as it acts as an efficient cooling agent in metal-poor gas (jllLPp . 
Indeed, a relatively large fraction of metal-poor stars below [Fe/H] = —2 have [C/Fe] > 1, as discovered by the 
HK and the Hambu rg/ESO surveys. The fraction of th e se car bon-enhanced metal-poor (CEMP) stars is estimated 



to be ~ 10 — 20% ( Cohen et al. . 2005 : Lucatello et al. . 2006f) . a fraction which tends to increase with decreasing 
metallicity. About 80% of the CEMP stars have been shown to be enhanced also in the s-process elements. This 
subgroup is accordingly named CEMP-s. There is now convincing evidence that the vast majority of the CEMP-s 
stars are members of binary systems and that they, like the more metal-rich CH and Ba stars, likely obtained their 
peculiar abundance pattern by mass transfer from an evolved binary companion that is now extinct (e.g., Fujimoto 
et al, 2000; Lucatello et al, 2005b). 

What about the remaining 20% of the CEMP stars? The majority of these stars are either enriched in r-process 
elements, but no s-process, (CEMP-r, i.e. CS 22892 — 052, Sneden et al, 2003) or they show no overabundance in 
n-capture elements altogether (CEMP-no). Even though attempts have been made to explain the abundance patterns 
of all CEMP stars with the mass transfer scen ario (e.g., Komiya et al, 2007; Masseron et al, 2010), it is not clear 
whether this is a necessary condition. Recentlv. ICooke et al ( 2011 ) reported on a metal-poor ([Fe/H] ~ —3) damped 
Ly-a system in the spectrum of the quasar J0035-0918, showing a significant enhancement in carbon. This finding 
suggests that extremely metal-poor star-forming gas could indeed be C enhanced, in which case true CEMP stars 
may exist without having to originate from a mass-transfer event in a binary system. Interestingly, three of the 
four ultra metal-poor stars (jChristlieb et al\ . [iooil : iFrebel et al\ . [ioost iNorris et al\ . [2007h found below [Fe/H] = -4.5 
are consistent with being CEMP-no stars. The only exception is the [Fe/H] = -4.7 star SDSS J102915 + 172927 



([CafFau et ali . 120111 ) which shows no evidence of CNO enhancements a nd is believed to have formed out of gas cooled 



by thermal dust emission ( Klessen et aZ.I . [2OI2I: ISchneider et al . 20121 ). Instead of treating these CEMP-no stars as 
individual, unique objects, it could, alternatively, be worthwhile to explore what one may learn if they were treated 
as 'normal' stars in the general context of chemical evolution. If the gas out of which the EMP stars were formed 
experienced a period of low or delayed star formation, e.g., due to negative feedback effect s from the first generations 
of stars, a small population of very C-enhanced, ultra metal-poor stars is to be expected (jKarlssoni 120061 ). Figure \W\ 
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TABLE III Redshift - lookback time. 

^ Ag6^ Lookback time Upper stellar mass^ 
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2.21 


11.58 
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1.018 
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5.98 


7.80 


0.934 


1.130 


0.434 


9.21 


4.57" 


1.089 


1.313 



" Age of the Universe at redshift of z. Calculated from the combined data (Maximum Lilselihood) of WMAP 7-year+BAO+-ffo 
(Komatsu et ai, 2011). 

Highest mass of still surviving stars that formed at redshift of z. These stars have a lifetime equal to the "lookback time". The first 
column gives the mass of stars with metallicity Z = 10~'^'^ Zq (which is taken as a proxy for Z = 0) while the second column gives 
the mass of stars with solar metallicity. The hyphens indicate that no stars existed at the time of the Big Bang. Data are taken from 
http: //stev.oapd inaf.it/cgi-bi n/cmd_2. 2 ( Marigo et al.„ ,2008,) . 
'^'Age of the Sun l|Bonanno et a i. 2002). 



shows the predicted distribution of EMP stars (grey and black shaded areas) in the [C/Fe] — [Fe/H] plane. Evidently, 
the predicted fraction of CEMP stars increases with decreasing [Fe/H], which mainly is due to the fact that low-mass 
star formation is assumed to be inhibited in regions with low contents of C (and O, see Bromm and Loeb, 2003; 
Frebel et al, 2007). In these models, the predicted fraction of CEMP-no stars below [Fe/H] = —2 is roughly 1 — 7%, 
depending on which set of stellar yields is used. This is consistent with the observed fraction of 2 — 5%. Note that 
this prediction would be only marginally affected by the presence of ultra-metal-poor stars formed by the dust-cooling 
channel (see Sec. IIII.DI) . unless this channel will be found to completely dominate over the fine structure line-cooling 
channel. As a result of the delayed star formation, the current model (Fig. \W\\ is also able to explain the apparent 
downturn of the number of stars in the metallicity range —5 ^ [Fe/H] < —4. The deficit of stars in this region should, 
in this scenario, rather be interpreted as an extension or a stretching of the Galactic halo MDF below [Fe/H] ~ —4. A 
small population of mega metal-poor stars ([Fc/H] < —6) is also predicted to exist as a result of an initial enrichment 
by electron capture SNe (Sec. IIII.B.2p . 

Chemical peculiarities which are a manifestation of mass transf er in binaries prove tha t it is difficult to match 
models of first star yields to the measured abundances of EMP stars (jMcWilliam et all [200i). This poses a potentially 
serious problem, as a non-negligible f raction of stars formed at all epochs may be formed in binary or multiple systems 
(|Stacv et aLl . l2010l : lTurk" et ali . liooi). Thus, the CEMP-s stars are not giving us a direct reading of the nucleosynthetic 
yields of the first stars. In particular, the light elements (e.g. CNO, n-capture elements, and possibly the a-elements) 
may simply reflect the yields of nuclear fusion in a massive companion (e.g. AGB), dredged up and transferred to 
the surviving companion. If the massive companion was even higher mass, the heavy elements of the survivor may 
reflect the nucleosynthetic yields at the time of the explosion rather than in a slow wind or Roche lobe overflow. This 
is a more complex situation than the more direct case when a star's abundances simply reflect the ISM chemistry 
at the time of collapse of the parent molecular cloud. In the latter case, the star's abundances are a true reading of 
the chemistry over some mass scale (i.e., the mixing mass or molecular cloud mass) in the early Univ erse. I n present 
day star clusters, chemical homogeneity to a high degree is clearly observed ( De Silva et aLl . l2007al lb. 20061) . We will 
explore this promising possibility of separating an original chemical signature from one that was imposed by a binary 
companion in more detail in Sec. IVIIII 



D. Constraining the primordial IMF 

1. The low-mass end of the IMF 



As the vast majority of the primordial stars presumably were massive objects with lifetimes 10 Myr, we do 



not expect m e tal-fr e e stars to exist today. This conclusion is well co r roborated by observations (iBeers et al. 



Carnev et all Il996t IChristlieb eTall . l2002l 12004 iFrebel et all l2005t iRvan and Norrisl llQQlt ISchorck et al 



1992; 



2009; 



Yannv et aLl . 120091 ). Due to a number of observational constraints, the search for metal-free stars has concentrated on 



the solar neighborhood and the Galactic halo, leaving the presently unreachable Galactic bulge unexplored. However, 
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the possibility that metal- free stars are exclusively hiding deep within the bulge seems slim (Brook et al., 2007, see 
also Sec. IV.B[) . Another concern has been raised that stars may sweep up enriched gas as they orbit around the 
Galaxy, thereby concealing their true nature, but this also appears unlikely (Frebel et al, 2009). To date, not a single 
metal-free star in known to exist in the Milky Way, or in any other galaxy of the Local Group that is observable with 
current facilities . This places direct constraints on the p rimordial IMF at the low- mass end (|Bond . 1981 : Karlsson , 
2005 : Oev . 20031 : [Salvadori et al. . 2007t [Tumlinsonl . |2006[) . and indicates that metal-free stars with masses < 0.8 Mq 



never formed (see Table HIT)) . 

CEMP stars may provide additional evidence for a higher mass scale of the IMF at early epochs. As already 
discussed, the dominant fraction of the CEMP stars are believed to originate from binary mass transfer. A high 
GEMP fraction like the one observed can be understood if the formation of single stars with masses ^ 0.8 A^© was 
suppressed at metallicities below [Fe/H]~ —2.5, which results in an IMF with a higher characteristic mass, mc, than 
that of the present-day one (e.g., Lucatello et al., 2005a; Komiya et al., 2007; Tumlinson, 2007b). For EMP stars, 
characteristic masses in the range 1 ^ rric/ M.q ^ 10 are found in the literature. This result is broadly consistent 
with early star formation being dependent on metallicity. It appears, ho wever, that metalli city may not be the 
sole parameter controlling the fraction of GEMP stars. Below [Fe/H]= —2, Frebel et al\ ( 20061 ) found an increasing 
GE MP fraction with i ncreasing distance from the Galactic plane in the Hamburg/ESO data. This spatial variation 
led Tumlinson (l2007bl) to argue for a time-dependent component of the IMF, coupled to the GMB (see Sec. IIILPp . 
The observed spatial variation of the fraction of GEMP(-s) stars can be explained in terms of a time-dependent IMF, 
where the Jeans mass is set by the decreasing temperature of the GMB, if the stars closer to the Galactic plane (i.e., 
thick disk stars) were formed later than the stars further out in the halo. Although plausible, the origin of the spatial 
variation of the GEMP star fraction merits further investigation. The fraction of binaries in a single stellar population 
could have been quite different at early epochs and may have varied with the star formation efficiency. If so, that 
would affect the fraction of GEMP stars. Nonetheless, irrespective of the precise explanation for the spatial variation 
of the GEMP fraction, the existence of such a variation can be utilized to obtain additional information on the origin 
of the Galactic halo and to help us discern the connection between the halo and the satellite dwarf galaxies. For 
example, a different GEMP fraction in dwarf galaxies as compared to that of, e.g., the inner halo, would be suggestive 
of different origins of the stars in the two type of systems. 



2. The high-mass end of the IMF 

As regards the high-mass end of the primordial IMF, theory and observations tend to disagree. To put constraints 
on the high-mass end, we must resort to measurements of the nucleosynthetic signature of the first generations of 
stars, preserved in the second and higher-order generations. 

As discussed in Sec. IVI.Al the emerging chemical signature of EMP stars in the Galactic halo indicates that the 
early stellar generations contained massive stars in a similar mass range and relative number fractions as the present- 
day IMF. There is no clear evidence of a significantly different shape of the IMF at the high-mass end. However, let 
us ask the question: assuming that primordial very massive stars were able to form and that at least some of them 
where able to explode as PISNe, in what metallicity regime should we expect to find their chemical signature? Due to 
their very short life-times, PISNe must have been the first stars to enrich the ISM in metals. In the classical picture 
of chemical evolution, an initial enrichment by primordial PISNe would generate a metallicity floor out of which the 
second generation, low-mass stars were able to form. Naturally, these stars would be the most metal-poor stars to 
be found in the Galaxy and would show a unique PISN signature, characterized by, but not limited to, a pronounced 
odd-even effect, highly super-solar [Si/0] and [S/G] ratios, and a lack of r- and s-process elements. Such a signature 
is, however, not observed, neither in the EMP stars, nor in the ultra metal-poor stars. These stars all show chemical 
signatures more resembling those of normal core collapse SNe (see, however. Sec. IVI.GI for a further discussion as 
regards the ultra metal-poor stars). This has been taken as an indication that very massive stars were exceedingly 
few in the early Universe (e.g., Tumlinson et al., 2004; Ballero et al., 2006), if not altogether absent, in contrast to 
the predictions. 

A different conclusion may, however, be drawn if the instantaneous mixing approximation assumed in homogeneous 
chemical evolution modelling is relaxed and the enrichment by SNe instead is allowed to occur locally, where the 
ejecta are mixed relatively slowly with their surroundings b y turbulent diffusion. In order to simulate the stochastic 
enrichment by primordial PISNe in a cosmological context. iKarlsson et al\ ( 20081 ) followed the initial star formation 



and chemical enrichment of collapsing atomic cooling halos. In the primordial star formation mode, only massive and 
very massive stars were able to form. As the gas became metal enriched, the primordial mode was suppressed and 
replaced by a normal Pop II star formation mode in which only low- and high mass stars were allowed to form. The 
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FIG. 17 Predicted distribution (shaded area) of low-mass stars predominantly (> 90%) enriched by PISNe in the [Fe/Ca] — 
[Ca/H] plane. Note that this is only a very small subpopulation of the full sample of stars in the simulation. The dashed, thick 
and thin lines indicate, respectively, the la (innermost 68.3%) a nd 3a (innermo st 99.7%) cosmic scatter of the full sample. 
Symbols denote various observations of Galactic halo stars. From iKarlsson et al\ (QOO^). 



larger explosion energies of the PISNe was also taken into account as they sweep up larger amounts of gas before 
merging with the ISM. 



The result is shown in Figure [TTl The predicted distribution of stars with a dominant contribution from primordial 
PISNe (gray-shaded area) is located around [Ca/H] ~ —2.5, which is significantly above the most metal-poor stars 
in the simulation (the distribution of the total population of simulated stars is indicated by the dashed lines) . These 
stars have, in fact, such high Ca abundances that a fair fraction of them may risk to remain undetected in surveys of 
metal-poor stars, such as the HK and Hamburg/ESO surveys, in which the metal-poor star candidates are ide ntified 
on the basis of the weakness (or absence) of the Ca ii K line. A similar result was found bv lGreif et all fooTdi) . who 
used three-dimensional cosmological simulations to study the assembly process of the first galaxies (cf. Salvadori et 
ai, 2007). For realistic estimates of the fraction of primordial very massive stars exploding as PISNe (say ^ 10% by 
mass in the primordial stellar population, Greif and Bromm, 2006), the fraction of low-mass stars bel ow [Ca/H] = —2 
with a dominant PISN signature is predicted to be very low, only ^ 3 x 10"'* ( Karlsson et ali . This may well 

explain the fact that such a star has not yet been found (see also Salvadori et ai, 200 7). Current observatio nal data 



are able to constrain the fraction of metal-free, very massive stars to ^ 40%, by mass (jKarlsson et ali . l2008l ) 



It should be pointed out that the newly discovered ultra-faint dwarf galaxies may be a particularly favorable site 
for searching for the PISN signature. If these galaxies originated from the atomic cooling halos they may have been 
enriched by, or even formed Pop HI stars (cf. Frebel and Bromm, 2010). Since a relatively large fraction of the 
member stars have metallicities below [Fe/H]= —2, the fraction of stars with a PISN-dominated signature may be 
larger with a factor of about 3, as compared to the corresponding fraction in the Galactic halo. 



We no te in passing that Eta Carinae and the Pistol Star are two examples of very massive stars in the Galaxy. 
Recently, ICrowther etai\ (|2010f ) also claimed to have identified sever al massive stars in exce ss of 150 A^0 in the young 
super star cluster R136 in the Large Magellanic Cloud. Furthermore, Gal- Yam et al. ( 20091 ) reported on an extremely 
powerful SN explosion which appears to have resulted from the death of a very massive star of about 200 A4q (see 
Smith et al., 2007 for another example of a very energetic SN explosion). 
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VII. CHEMICAL SIGNATURES FROM LOW MASS GALAXIES 
A. Dwarf galaxies 

Low-mass galaxies in the Local Group have come increasingly into focus as ideal laboratories to study the early 
enrichment history of the universe. Here, we first introduce the relevant phenomenology, and proceed to discuss what 
constitutes a bona fide 'first galaxy'. These systems, in turn, could serve as the host for the formation of the first 
low-mass. Pop II stars, thus allowing us to make the connection with the nucleosynthesis provided by the first stars. 



1. Simplified chemical laboratories 



In recent years, we have come to realize that the Galactic halo is an amalgam of many systems and fragments 



(|Helmi l2008t iTolstov et all . l2009l ) 



Rather than concentrating on the general halo population, we anticipate tha t 
future studies will target metal-poor stars in halo subsystem s. These includ q^^ | dwarf spheroidals (Mateo . 



ultra-faint dwarfs (IKirbv et al. 



stellar associations ( Walsh et al 



2008t [Simon and Gehal . l2007t) . stellar s treams Ichou et all [2010t llbata et al 
2007 ). and satellites to dwarf galaxies ( Belokurov et al . 20091 : Coleman et al 




So to what extent can we learn about the first stars from these systems? 

Initially, things did not look good for fossil signatures in dwarf galaxies. Se veral groups pointed out that stars 
enriched by the f irst stars are likely to be found at the centres of galaxies today ( Bland- Hawthorn and PeeblesL 2006t 
iTuml msonl . I2OIOI) although other s have shown that first stars can be delayed to z ^ 3 such that prese nt-day dwarfs may 
show the signatures of first stars ( Brook et al. , 2007 ; Scannapieco et al. , 2006t iTornatore et al. , 2007 ) . Observationally, 
the classical dwarfs Sculptor, Forna x, Carina and Sext ans were compared to the halo stars, but no evidence for 
stars below [Fe/H]= —3 was found (jHelmi et all . I2OO6'). But a re-calibration of the Ca 11 triplet (C aT) lines has 
revealed numerous metal poor stars below [Fe/H]= — 3 although many of these have yet to be confirmed (jKirbv et all . 
2008 : Norris et all . l2008l Starkenburg et al . 2010[ l. Seemingly like stars in the outer Galactic halo (e.g. Roederer , 



2009), evidence is now emerging t hat star-to-star abundance variations exist in dwarf galaxies ( Feltzing et all . I2OO9I: 
Fulbright et~ai\ . |2004 iKoch et all I2OO8. ). We now explore the prospect that dwarf galaxies are important sites to 



learn about the chemistry of the first stars after all. 

The least massive dwarf gala xies, if these can be unam biguously identified, may provide the cleanest signatures 
of the yields of the first stars ( Frebel and Bromm . 2010f ). The very low rate of star formation measured in the 
faintest dwarf galaxies presumably indicates that there have been relatively few enrichments by supernovae over the 
lifetime of the dwarf. This means that the mean metallicity is expected to be lower thereby giving a higher fraction 
of metal poor stars. Furthermore, we may expect that the abundance scatter is large at low metallicity consistent 
with the small number of enrichments for the low baryon mass fraction. Complicating factors are the effects of 
feedback and environment. The well known trend of decreasing metallicity with declining galaxy luminosity has been 
argued to be a consequence of metals lost through galactic winds (see Dalcanton, 2007), but others have argued that 
this largely reflects the lower star formation efficiencies in smaller galaxies (see Tassis et al, 2008). Presently, it is 
unclear what impact this will have on the overall abundance scatter. Furthermore, the environmental influence of 
dynamical stripping is also uncertain. If faint dwarfs were the stripped- down cores of larg er dwarfs, would these 
look chemically distinct from present-day dwarfs of the same luminosity? iKirbv et all (12008^ demonstrate that the 
luminosity- metallicity relation is present down to log{L/CQ) ~ 4 but it is unknown wheth er this reflects a simp le 
mass- metallicity dependence or a more complex process involving, e.g., the baryon fraction (|Salvadori et al\ . \201(xi . 



2. Minimum galaxy mass 



What then is the minimum mass a galaxy can have (Bromm and Yoshida, 201 1)? This question has been asked 
numerous times over the years in different conte xts (see, e.g., Tols t oy et al, 2009) and i s parti cularly relevant today 
with the discovery of ultra- faint dwarf galaxies ( Simon and Geha . 2007; Willman et al . [20051 ). In part, the answer 
depends on the epoch of galaxy formation. If dark matt er is made up of weakly interacting massiv e particles (mw ~ 
100 GeV), it can fragment down to Jupiter mass scales ( Bertschingei . 20061 : Diemand et al . 2005 ). But in order for 
the IGM to overcome thermal pressure and accrete onto dark matter, the mass of the minihalo must have exceeded 



We have not included globular clusters in this list. W hile they con stitute a halo population, they have remarkably high levels of 
enrichment for systems that are believed to be ancient IGratton et all . [2004) . 
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10^ A4q before reionization (|Loeb et all l2008l) . During the reionization epoch, in order for the local gas density 
to survive heating through photoionization ( which suppresses star formation), the minimum halo mass exceeded 

These simple calculations based on ionization balance 



10^ Mq (|Haiman et g/J . Il996t llkeuchi Il986l: [Recs 



are supported by 3 D hydrodynamical simulations (jGnediil 120001 : Navarro and Steinmet j . 1997 ; Okamoto et al. , 20081: 



Quinn et a/.l . ll996l ). In contrast, others have argued that star formation was diminished but not ful ly suppressed, such 
that galaxies with masses below 10^ A^f?) form stars dur i ng an d immedi a tely a fter reionization ([Bovill and Ricotti . 



2009t Gnedin and Kravtsov . 2006t Salvadori and Ferrara , 2009() . Ricotti ( 2009t ) has argued for late-phase accretion 
in dwarfs as recently as z = 1 — 2. Independen t evidence of this may be the starbursts spaced by a few Gyr seen 

A I2OO4I : pRlstov et all . l2009[ ) : these bursts were 



in the star formation histories of dwarf galaxies (jGrebel and Gallagheii 
presumably triggered by accretion or re-accretion of cold gas, although a perigalactic passage is not ruled out. 

Assuming that gas accretion can occur, we will take as the working definition of a "minimum mass galaxy" a dark 
matter halo whose baryons survive the impact of (at least) one supernova explosion. After wards, the enriched gas 
should be able to cool, collapse and form (at least) one subsequent generation of stars. Recently, iBland -Hawthorn et al\ 
(|2010cl ) demonstrated that a clumpy medium is much less susceptible to SN sweeping (particularly if it is off-centred) 
because the coupling efficiency of the explosive en ergy is much lower than for a diffuse interstellar medium. With the 
aid of the 3D hydro code Fyris ( Sutherland . 2010t ). they show that baryons are retained and stars are formed in dark 
matter haloes down to 3 x 10^ A^©- These systems are expected to have distinct chemical signatures which may be 
detectable in the near future, either from direct observations of dwarf galaxies, or from quasi-stellar object (QSO) 
sight-line spectra due to the presence of dwarfs around distant galaxies. 

This limiting case opens up a profitable line of enquiry , part icularly when we consider the chemical signatures 



produced by the discrete supernova events ( Karlsson et al . 2010f) . While the minimum mass galaxies may not con- 



tribute a large fraction of baryons and dark matter to the formation of galaxies, they are likely to carry important 
information on the stars that were formed at the earliest cosmic times. The possible connection between the min- 
imum mass galaxies and the faintest dwarfs in the Local Group must be further investigated. In fact, given that 
stars below [Fe/H 



Starkenburg et al 



3 a re detected in the dwarf galaxies (jFrebel et al. . l2010btlKiri)v et a/.l . [iool iNorris all . l20ld 
we could, without fully realising it, already be looking at the very formation sites of the 



first stars (cf. White and Springel, 2000). 



B. Damped Lyman-a systems 



An interestin g deve l opment is the rece nt discovery of very metal-poor DLA systems along QSO sight lines 
( Penprase et al . 2010l : Pettini et al. . 2008 ). If these are protogalactic structures that have recently formed from 
the IGM, there is the prospect of identifying the chemical imprint of early generations of stars (e.g. Pettini et al, 
2002). Interestingly, two of these sy stems with [Fe/H] ^ —3 may bear the hallmarks of early stellar enrichment 



(|Cooke e^a/lbOllI: 



Erni et al. 



200&) . 



First, lErni et all (|2006[ ) identify a DLA toward the QSO Q0913-I-072 (z = 2.785) with an iron abundance charac- 
teristic of the IGM at that redshift. The C, N, O, Al, Si abundances show an odd-even effect reminiscent of the most 
metal-poor stars in the Galactic halo. This pattern is created in models where the neutron flux is low (e.g. Heger and 
presumably due to the low overall metal abundance. A more striking signature is the strong [N/H] 



depletion which 


Pettini et al. 


IGM. Erni et al. 


(2006) argue 



Pettini et ali (|2002l ) has argued is further evidence for a system that has recently formed from the 



(|2006l ) argue that the abundances are in good agreement with 10—50 M© zero metallicity Pop HI 



models. 

Secondlv. ICooke et all (|201ll ) identified a DLA towards the QSO J0035-0918 (z = 2.340) with a remarkably strong 
[C/Fe] enhancement and with the same odd-even effect in light elements. This result is particularly striking because it 
is rem iniscent of the recent discovery of CEMP stars in the Galactic halo ([Beers and Christliebl . [2005t iLucatello et all 
I2OO6I ). They infer the total mass of neutral gas to be ^ 3 x 10^ Mq within a linear scale of ^ 100 pc. The inferred 
amount of carbon is consistent with one or two SN enrichment events from primordial, or very metal-poor high mass 
stars. At face value, these observations suggest that at least some fraction of the CEMP stars, in particular those 
which are not s-process enhanced (see Sec. IVI.CI). may have formed out of carbon-enhanced gas, and are therefore, 
in a sense, truly C-enhanced stars ( KarlssonT 2006[) . 

The importance of such observations is two-fold: (i) they provide a chemical signature on a mass scale that is 
orders of magnitude larger than the mass of a stellar envelope, and is therefore expected to be more representative 
of early star formation; (ii) they allow us to probe in situ enrichment at early cosmic epochs. The fact that we 
recover chemical signatures that are strikingly similar to those ob served in ancient halo stars is a strikin g example of 
the powerful synergies between near-field and far-field cosmology ([Freeman and Bland-Hawthorni 120021) . In the next 



38 



section, we look at another example of a chemical signature on comparable mass scales. 



VIII. THE IMPACT OF STAR CLUSTERS 



As we have seen, the chemical infor mation arising from the most metal-poor stars is very difficult to unravel 



(jKirby et all 120081 : iNomoto et al\ . |2005[ ). One possibility is that not all of the stars are providing us with an unam- 



biguous 'reading' of the early enrichment of the primordial interstellar gas. As mentioned in Sec. IV I ^C j a sign ificant 
fraction of extremely me tal poor stars appear to have undergone mass transfer with a companion ( Lucatello et al\ . 
l2005bl : lRvan" et all I2OO5I) which complicates any attempt at inferring the progenitor yields for elements such as CNO 
and light a elements. Binarity m ay explain why the elemental abundance s of th e most metal poor stars defy a clear 



explanation at the present time (jJoggerst et all l2010at iMc William et all [200i). It is therefore imperative that we 



can distinguish the effects of binarity from the effects of inhomogeneous mixing, where we may stand a better chance 
of learning about the yields of the fi rst sta. rs. 

Recently, ISland-Hawthorn et al\ ( 2010a ) demonstrate how this issue may be resolved by incorporating a missing 
ingredient into existing models of stochastic chemical evolution - the formation of star clusters. If the surface 
abundances of a star are altered, e.g. by a mass transfer from an evolved binary companion, they no longer constitute 
a bona fide imprint of the chemical state of the gas cloud out of which the star once was formed. However, if one 
could have access to a large number of stars formed out of the same gas cloud and with identical surface abundances 
to begin with, one could afford a fraction of these stars to be 'polluted', without suffer the loss of information on the 
parent gas cloud. This is where the star clusters come in. If clusters were able to form in the early Universe, they 
should be incorporated in the models to allow for a correct interpretation of the observational data. 



A. The distribution of cluster masses and chemical homogeneity 

In the present-day Universe, most stars are born in a single burst within compact clusters and stellar fragments, 
rather than in isolation. This fact is well established in the local Universe (Lada and Lada, 2003) and a growing 



body of evidence suggests that it is likely to be true also at the time of the first galaxies ( Clark et al . 1200811201131 : 



Karlsson efah . I2OI2I : IStacv et al\ . bOld iTurk et al\ . l2009h 



In order to derive the impact of star clusters on the abundance plane, we must consider the progenitor mass 
distribution of star clusters. It is now well established (e.g. Larsen, 2009) that star clusters have a range of masses 
extending from a minimum mass (Mmin) to a maximum mass (Mmax)- Although the particular situation for dwarf 
galaxies is not clear, we know from the cloud mass distribution that the range of masses covers many decades (e.g. 
Escala and Larson, 2008). The birth distribution of stellar clusters is known as the initial cluster mass function 
(ICMF) and is assumed to have the form 



AN/dKU = x(M,) = xoM; 



(17) 



where M* denotes the cluster mass 
slope of 7 w 2 in most environments 



At least in the nearby Universe, t he observations may support a universal 

equal mass per logarithmic bin ([Elmegreen . I2OIOI: iFall et all J2005, 2oor 



Lada and Lada . 20031) . although a flatter slope may be more applicable in the early Universe ( Bland-Hawthorn et all . 



2010a[ ). However, the slope cannot be too flat. Since the star-to-star scatters, particularly for the r-process elements, 
decrease with decreasing 7, the slope must be steep enough to ensure consistency with the observed scatters (e.g., 
F rancois et al, 2007). 

De Silva et al\ ( 2007al [bl. l2006h have shown that both old (~ 10 Gyr) an d intermediate-age (~ 1 Gyr) op en clusters 
are chemically homogeneous to a high degree (A[Fe/H] ^ 0.03 dex). iBland-Hawthorn et al. ( 2010bl ) provide a 
condition that explains why chemical homogeneity is expected in most star- forming clusters. For a star cluster to 
be chemically uniform, it must form from a gas cloud that has collapsed and finished forming stars before the first 
supernova goes off within the cloud. A 'homogeneity condition' is given that depends only no the column density 
through the cloud for a given cloud mass. A typical open cluster is expected to be chemically uniform up to about 
10^~^ M0, and a globular cluster is expected to be uniform up to 10^"'^ Mq. The homogeneity condition has an 
important consequence for the distribution of stars in the abundance plane, and the effects should be observable even 
in the limit of only a few data points. 



Apart from a few light elements, the same holds true for globular clusters llGratton et aZ.I , |2004) . although few systems (e.g., a;Cen and 
NGC 1851) recently have been found to show evidence for more than one burst of star formation (Lee et al., 2009). 
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As a result of stars being formed in chemically homogeneous clusters, "clumping" of stars in abundance space, say, in 
the [Eu/Fe] - [Fe/H] plane, is to be expected. It is important to understand that the majority of the clusters considered 
here are supposed to have been dissolved into the diffuse field. However, while the stellar kinematical information 
is lost, the chemical information is still preserved and the prospects of reconstructing clusters in abundance space 
are high, as long as the elemental abundance determinations are accurate enough (i.e., A [Fe/H] ^0.1 dex). This is 
the key point. If detected, the most metal poor abundance groupings are likely to reflect the conditions in the gas 
at the onset of star formation. The spectra of these stars can be added without loss of information to produce a 
more accurate measurement of the progenitor conditions. This ab undance measurement is averag ed over a substantial 
amount of gas and is therefore not subject to mixing anomalies (jKarlsson and Gustafssonl . 120011 ) or mass transfer in 
binaries (e.g. Lucatello et ai, 2005b). 

A clean signature of clustering, particularly at low metallicity, is important because it indicates the presence of 
massive star clusters in the early Universe and conceivably provides a constraint on the mass of the first systems. 
Particularly strong clumping in abundance space, for a given number of stars in the sample, would indicate a highly 
flattened ICMF, or high mass cut-offs (Minin,Mniax). This could herald the onset of the formation of massive star 
clusters in dwarf galaxies (e.g. Bromm and Clarke, 2002). If the star formation efficiencies were low at that time, this 
may require supermassive gas clouds (^ 10^ -Mq) to have formed even at the earliest times (I Abel et al . 2000), possibly 
consi s tent with the regular o c currence of massive star-f orming clumps at high redshift (|Elmegreen and Elmegreen . 



2005 : Forster Schreiber et al . 20061 iGenzel et al. . 2006 ). Conversely, if such clustering was not observed, then we 



would infer that the slope of the early ICMF is steep, or the maximum cluster size is relatively small compared to the 
present day. But the observed scatter in the abundance plane would need to be consistent with the non-detection of 
clustering. 



B. The Galactic halo 



Clumping in abundance space can also be used as an independent probe of galaxy formation. The detection of 
clumping in the Galactic halo would give us unique insight into how the different components of the Galactic halo 
were formed and whether the building blocks of the halo possibly differed, in terms of how they formed stars, from 
the dwarf galaxies that presently orbit around the Milky Way. However, since the stellar mass of the halo is large 
as compared to the mass of a typical cluster, the number of clusters building up the halo is relatively big (cf. Sec. 
IVIII.C|) and clumping may be difficult to detect. 

Let us illustrate with a specific example. The large SDSS/SEGUE surveys (see Sec. IV. Al) contain data on hundreds 
of thousands of stars, most of which belong to the Galactic thick disk and halo. Although the spectral information 
is limited, these surveys have unveiled a new multi-component Galaxy rich in substructure and whose formation 
and subsequent evo l ution generally is con s istent with a hierarchical build-up from smaller systems ( Bell et al. . 2008t 



Carollo et al. 



( Yannv et 



20ld llvezic et all l20Q8[ ). Among the SEGUE stars, - 60,000 are G dwarfs 



. l201fll : Id^ 

20091 ). Only a minor fraction of these stars likely belong to the halo. Assuming a three-component, 
smooth model of the Galaxy (see, e.g., de Jong et al., 2010), we estimate that about 35% of the SEGUE G dwarfs 
are halo stars, which corresponds roughly to 21,000 stars. 

Given that they can be identified unambiguously, how much of the total (stellar) halo mass is probed by these 
G dwarfs? If the stars are completely phase mixed, they will probe the entire initial formation mass Mhaio.form — 
1 X 10^ A^0 of the Galactic halo. An extreme lower limit would be to assume that the stars only probe the instant 
volume of the survey. We can, however, do slightly better by assuming that all stars follow circular orbits around the 
Galaxy and integrate the mass over the spherical shell whose boundaries are defined by by the inner- and outermost 
G dwarf orbits. Again, making use of a smooth model of the stellar halo, and assuming that the G dwarfs are found 
within 6 — 14 kpc from the Galactic center, the mass fraction probed by the SEGUE halo stars is estimated to ^ 0.1. 
Since halo stars usually have non-zero eccentricities, this estimate should be a lower limit. 

The two scenarios are displayed in Fig. [18] (top and middle panels). The left- hand-side panels show the distribution 
of simulated Galactic halo st ars in the [Ca/Fe1 - [Fe/H] pl ane. These are generated by the stochastic chemical 
evolution model presented in iBland-Hawthorn et al. ( 2010al ). assuming a present-day-like ICMF with 7 = 2 and 
cluster masses between 50 and 2 x 10^ Mq- The corresponding right-hand-side panels show the excess amount of 
detected clumping (blue lines) as corn pared to spurious, random clumping (black lines). The group-finding algorithm 
EnLink (jSharma and Johnston) . l2009h is used to produce the clumping statistics. EnLink calculates the number of 
clumps having a given number of stars within them, both for the wanted distribution and for a "smooth" background 
distribution. A comparison between these two numbers indicate whether excess clumping is present in the wanted 
distribution. The "observational" uncertainties are set to 0.01 dex both in [Ca/Fe] and [Fe/H]. Although beyond 



40 



1.2 



0.8 

^ 0.4 
O 



0.0 



• 














,100 



10 A, 

CD 




[Fe/H] 



FIG. 18 The possibility of detecting star clusters in the Milky Way halo. The left panels show the predicted distribution of 
21,000 stars in the [Ca/Fe] - [Fe/H] plane. The intrinsic "observational" uncertainties are set to 0.01 dex in both axes. The 
right panels indicate the corresponding degree of excess clumping (blue lines) as compared to spurious clumping in an "infinitely 
massive" model where not two stars originate from the same cluster (black lines). T he clumping statistics is calculated for 
stars below [Fe/H] — —2.5 (white areas) using EnLink (|Sharma and Johnstonl 120091 '). The parameter G{> S) is defined as 
the number of clumps with > 5* stars relative to the background (S is defined as a relative overdensity) . The stars in the top 
panel probe the entire initial stellar mass of the Galactic halo (Mhaio.form = 1 x 10^ Mq) while the stars in the middle panel 
only probe 10% of the halo mass (i.e., 10% of the number of star clusters), which is comparable to the mass probed by the 
halo G dwarfs in the SEGUE survey. In contrast to the top panels, clumping is clearly observed. In the bottom panel, the 
stars from the middle panel are additionally impaired by an observational uncertainty of 0.2 dex in both abundance ratios. As 
a consequence, the clumping is almost entirely washed out. The red dots in the left panels denote stars strictly enriched by 
material from a single SN. These stars may not necessarily be the most metal-poor stars in the sample. 



reach for present-day abundance analyses, we choose such a small uncertainty in order to disentangle the observational 
"smearing" from the "smearing" due to under-sampling. Clearly, clumping is intrinsically difficult to detect for the 
larger halo mass Mhaio.form — 1 x 10^ -Mq, and for the given number of stars in the survey (top panels). In a lOx 
smaller halo mass, however, excess clumping is present (middle panels). Note that the clumping of stars is visible for 
the naked eye in the middle left panel. The clumping statistics displayed in the middle right panel clearly shows the 
presence of excess clumping. The bottom panels show the resulting amount of excess clumping taking into account 
the observational uncertainty of the SEGUE survey. A fiducial uncertainty of 0.2 dex in both axes will almost entirely 
wash out any excess clumping. This suggests that evidence of clustering will be difficult to detect in the SEGUE 
sample, even if stars did form in clusters also at the epoch of formation of the Galactic halo stars. 

We note in passing that stars strictly enriched by a single SN (red dots in Fig. [TE\i constitute a fraction of 1.5 x 10~* 
of the total number of halo stars in this particular model. The corresponding fraction below [Fe/H]= —2.5 is 1.3xl0~^. 
These fractions depend on several factors such as the IMF, the ICMF and the amount of turbulent mixing of the 
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FIG. 19 The expected amount of excess clumping (blue lines) in two different model halo samples. The left panel display the 
amount of clumping in a sample 10 times the size of the SEGUE halo G dwarf sample (~ 210, 000 stars, 0.2 dex uncertainty 
in both axes) while the right panel display the amount of clumping in a SEGUE-sized sample (~ 21,000 stars) for which the 
observational uncertainty is reduced to 0.05 dex in both axes. Excess clumping is detected in this sample, while the sample 
displayed in the left panel only show very little evidence of excess clumping. See caption of Fig. [18] for a definition of S and 
G(> S). 

ISM. These stars generally appear in a range of metallicities from —4 < [Fe/H] < —2.5 and may not necessarily be 
the most metal-poor stars in the sample. Note that these stars are enriched by normal core collapse SNe and are 
distinct from the 'second generation stars' predominantly enriched by the PISNe, which may show up at somewhat 
higher metallicities, as discussed in Sec. IVI.DI The identification of 'single-SN-star-clusters' will give us direct insight 
into primordial stellar yields. 

Figure [in] shows the amount of excess clumping predicted for a sample 10 times larger than the SEGUE halo G dwarf 
sample (left) and the amount predicted for a SEGUE-sized sample for which the observational uncertainty instead is 
reduced to 0.05 dex in both axes (right). Enhanced excess clumping is detected in the latter case. This shows the 
necessity of acquiring as small observational errors as possible, sometimes in favor of a larger number of stars, since 
crucial information may be buried in the details. If the abundance errors are worse t han about 0.1 dex, to improve 
our c hances of detection, we should look to use as many distinct elements as possible fiBland- Hawthorn and Freeman . 
20041) . The probability of detecting clumping will significantly increase with increasing number of elements, in par- 
ticular for n-capture elements like Eu and Ba, which are expected to show a large scatter in their abundance ratios 
(Fig. [20l see also discussion in Bland-Hawthorn et ai, 2010a). The number of stars must, however, be large enough 
to ensure an adequate sampling rate per cluster. 

In the (worst) case of complete phase mixing of the Galactic halo stars (i.e., a, by the survey, sampled mass of 
Msurvcy = Mhaio.form = 1 X 10^ ■Mq), wc uccd roughly 2.4 X 10^ stars to obtain an average sampling rate of one 
star per cluster (note that the most massive clusters will be sampled by 3> 1 stars while the least massive ones will 
generally fall below a sampling rate of one star per cluster), assuming a present-day- like ICMF with cluster masses 
in the range 50 — 2 x 10^ Mq and a slope of 7 = 2. This is 10 times more stars than the total number of stars in 
the SEGUE survey. Hence, depending on parameters like the ICMF and the amount of phase mixing, there is an 
optimal survey size for identifying evidence of clustering. However, in order to probe, e.g., 'single-SN-star-clusters', 
the low-mass end of the ICMF, or the enrichment by rare types of SNe, bigger samples may be required. 

C. Dwarf galaxies 

So far, we have mainly focused on the Galactic halo when discussing the possibility of finding evidence of ancient 
star clusters by searching for clumps or aggregates in the chemical abundance diagrams. What about the stars in 
dwarf galaxies? As suggested from the results in the last section (see Fig. [T8]) . the amount of excess clumping is 
expected to be larger for less massive systems, i.e., for a fixed number of stars in the sample, and the crowding should 
be less severe. Unfortunately, the larger distances to the dwarf galaxies make it a challenge to observe many stars. 

To emphasize the impact of clustering with low number statistics, i.e., we simulate the abundance measurements 
for a dwarf galaxy at a distance of 30 kpc (see Fig. [20]), as observed on 8m class and 30m class telescopes, respectively. 
We consider a galaxy with a stellar mass of 3 x 10^ -Mq, typical of a faint dwarf galaxy. This object has about 10^ 
stars, a luminosity of 10^ Cq and an absolute V mag of My = —7.6 assuming a Salpeter IMF. For echelle observations 
on an 8m telescope (Fig. [20]) . we assume measurement errors of 0.1 dex in both [r/Fe] and [Fe/H]. Since the mass 
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FIG. 20 A simulation of a targeted study of a nearby (D = 30 kpc) dwarf galaxy (M* = 3 x 10® -Mq) on an 8m (top) and 30m 
(bottom) telescope, assuming a limiting magnitude of V = 18 and V = 22, respectively. Stars are assumed to form in clusters 
according to an ICMF with a slope of 7 = 1.5 (left). The panels on the right show a corresponding galaxy with no clustering, 
for comparison. Simulated errors are 0.1 dex in both [r/Fe], r standing for any r-process element, and [Fe/H] in the top figures 
and 0.05 dex in the bottom figures. There is evidence of clustering at [Fe/H] < —3.0 from a sample of ~ 10 stars on an 8m class 
telescope; the clustering is easily detected in the 30m telescope experiment. Above [Fe/H]= —3 (gray sh aded area), c rowding 
becomes increasingly severe which makes clumping increasingly more difficult to detect. From iBland-Hawthorn et al\ l(2010a ). 



distribution of star-forming clouds in dwarf galaxies is not well established, we treat the slope of the ICMF in Eq. 
(|17p as a free parameter. For 7 < 2, the effects of clustering start to become evident. For the same experiment on a 
30m telescope (Fig. I20[) . we assume a general improvement in the atmospheric models and the experimental errors, 
and therefore adopt errors of 0.05 dex. The effects of clustering, which are easily seen, remain clearly visible even 
after a twofold increase in the measurement errors in both axes. 

In summary, the simulation in Fig.[20]is a powerful statement of the importance of multi-object echelle spectrographs 
on ELTs. It should be relatively straightforward to detect the decreasing abundance scatter and the imprint of early 
star clusters for individual dwarf galaxies. And, by specifically targeting subgiants, for which accurate ages can be 
determined, we will have an additional, independent measure of cluster membership. The reconstructed cluster must 
be a coeval population by definition. 
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1. Signatures of ancient star clusters: the case of the MDF 



Signatures of relic clusters does not have to be identified in multi-dimensional abundance space. Already in one 
dimension, "bump s" and irregula, r ities i n the MDF, particularly in the metal-poor tail, could be a indication of 
disrupted clusters. iKarlsson et al. ( 2012 ) report on a tentative detection of a relic star cluster in the Sextans dSph . 
They re-analysed the excellent abundance data of six very metal-poor stars below [Fe/H] = —2.5 bv lAoki et al. ( 20091 ) 
and found that at least three stars clumped together in chemical abundance space, particularly in Mg, Ti, Cr, Ba, and 
Fe. While two of the "clump" stars have similar [Na/Fe], the third star have a significantly lower [Na/Fe]. This could 
possibly be a manifestation of the Na-0 anti-correlation observed in Galactic globular clusters, and may suggest that 
these stars were once members of a very metal-poor globular cluster. The weighted iron abundance of the "clump" 
stars is [Fe/H] = —2.7, which would ma ke the parent cluster the most metal-poor cluster to date. Interestingly, in 
the re-calibrated CaT data of Sextans bv IStarkenburg et al. ( 2010[ ). there is a clea r sign of a "bump " in the MDF at 
[Ca/H ] ~ —2.65, very close to the mean Ca abundance of the three stars observed bv lAoki et all ( 20091 ). iKarlsson et al. 
(I2OI2I) made use of this "bump" to estimate the initial mass of the cluster to M^^init = 1.91q g(1.6lQ g) x 10^ Mq, 
assuming a Salpeter (or Kroupa, Kroupa, 2001) IMF. 

There is also evidence for a kinematical substructure, possibly a remnant star cluster, in the central region of Sextans 
(Kleyna et al., 20 04; Battaglia et al., 2 011; see, however. Walker et al., 2006). This substructure has a metallicity 



of [Fe/H] ^ —2.6 (jBattaglia et aZ.1 . 120111 ). which is close to the metallicity of the "clump" stars and a common origin 



would not be implausible. This lends further support to the hypothesis of the existence of a now disrupted star cluster 
in Sextans, whose signature can be detected both in velocity space, multi-dimensional abundance space, and in the 
MDF. 

Met al-poor environmen t of dwarf galaxies like dSphs provides an ideal base for searching for ancient clusters. In the 
data of lKirbv et al\ (|201l[ ). several MDFs exhibit "bumpy" tails, such as Leo II, Draco and Ursa Minor. Interestingly, 



the U rsa Minor dSph is known to possess a second concentration of stars roughly 14' from the center. iKlevna et al\ 
(|2003l ) claimed that these stars belong to a cold substructure, consistent with the remnants of a disrupted star cluster. 
It would als o be i nteresting to compare ensembles of MDFs for different types of gal axies. In a first at t empt . 



Karlsson et al. ( 2012 ) compared the MDFs of Sextans, Sculpt or, Fornax, and Carina d Sph ( Starkenburg et al. . 2010l ) 



with the MDFs of eight ultra- faint dwarf galaxies observed bv lKirbv et al. ( 2011 , 2008 ). They found an excess scatter 
between the MFDs of the dSphs, consistent with the presence of clusters distributed in mass according to the present- 
day ICMF. No such excess scatter was found between the MDFs of the eight ultra-faints. Together with an observed 
offset between the average MDF of the ultra-faints and that of the dSphs, the distinction in the degree of "MDF 
bumpiness" may suggest that the two groups of galaxies were formed in different environments. Again, this shows 
that the identification and reconstruction of disrupted star clusters in chemical abundance space is important, not 
only to gain a deeper understanding of star formation and stellar yields, but also of galaxy formation and the build-up 
of galactic stellar haloes. Any "bump" in the MDF of a dwarf galaxy need, however, to be confirmed/refuted as a 
cluster signature, e.g. in multi-dimensional abundance space. 



IX. PEERING INTO THE FUTURE 



It is an extraordinary fact that we can probe back to the earliest generations of stars from observations of the local 
Universe — this is what we mean by 'near field cosmology.' The fi rst stars were unique to their time: they were 



responsible for the first chemical e lements (IRyan- Weber et ali . 120061 ) and for reionizing the neutral fog of hydrogen 



that permeated the early Universe ( Fan et all 2002 ). Our goal in describing pre-galactic enrichment is to understand 
the formation of the first stellar generations and their yields. Did the first stars form in isolation or in groups? Were 
relatively few stars responsible for reionization or was it triggered by the collective effect of massive star clusters? Just 
what are the processes that govern star formation at extremely low metallicity? Is this exclusively the domain of the 
most massive stars, or can substantial intermediate and low mass stars form? In other words, did stellar populations 
observable today exist before reionization? We have few if any answers at the present time. While this is a field that 
is very much in its infancy, we believe there has been sufficient progress to warrant our review. 

The chemical signatures imprinted in old, metal-poor stellar populations make up a vast repository of information 
that takes us back to the first few hundred million years of the Universe. It is an invaluable tool that allows us to 
witness the very process of primordial and early star formation in unparalleled detail. But not only that, chemical 
signatures on galactic scales help us to understand the connection between local dwarf galaxies and those systems 
which once built up the Milky Way. A more coherent and complete picture of the galaxy formation process is now 
emerging from observations of nearby stars. 
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We anticipate a great deal of progress over the next decade, particularly in an era of extremely large, ground-based 
telescopes (25-40m) and a new generation of space telescopes. We anticipate that many of the great advances will 
come from more extensive observations of the near field, in particular, metal poor stars in the vicinity of the Galaxy. 
But we should not rule ou t the extraordinary . For example, at least one gamma-ray burst has been visible to the 
naked eye in recent years ( Bloom et ai . 2009() allowing for high resolution spectroscopy which can in turn provide 
important constraints on the local stellar populations (e.g. Castro-Tirado et ai, 2010). High redshift quasars beyond 
2; ~ 3, even with the added bo ost from gravitational lensing, are likely to be too faint to provide useful absorption line 
spectra of intervening galaxies ( Glikman et al. . 2008[ ). We may obtain future constraints from rest-frame UV ionization 
diagnostics of high-redshift sources. Several aut hors have identified a population of CNO-rich galaxies being iii ternally 
ionized by star clusters with a top-heavy IMF (jFosbury et al. . l2003HGlikman et ad liool iRaiter et aLl.l2010l). 

Th e study of the high redshift Universe is complementary to studies of the local Universe (jBland-Hawthorn and Freeman 
20001 ). For example, measurements of the integrated spectral signature of Pop III stars is an exciting prospect although 
separating the signal fro m the infrared foregroun d is difRcult (see, e.g., Raue et al., 2009). Even though tantalizing 
progress has been made (jKashlinskv et aLl . l2007[ ). we must await the advent of next generation telescopes, such as 
the JWST and ALMA to be able to resolve the first galaxies, and detect individual sources and their birth sites 
(jWalter and Carillil 120081 ). It should be stressed that the level of detail on the first stars and galaxies, buried in the 
chemical abundance patterns of old stellar populations, can never be reached by studying the direct light from these 
objects themselves, nor by studying the high- 2: IGM. The information needs to be extracted from resolved stars. The 
modeling of the build-up of chemical elements over cosmic time, in particularly at early stages, is absolutely essential 
to the interpretation of this information. 

But our expectation is that th e most important insig hts on pre-galactic enrichment will come from the near field. 
We may be in for some surprises. iKarlsson et al\ pOOST ) has described how a significant fraction of second generation 
stars that formed out gas enriched by PISNe are, contrary to what may be expected, predicted to be found at 
relatively high metallicities around [Fe/H] ~ —2.5. For this reason, they may not have been picked up in existing 
surveys of metal-poor stars, like the HK and Hamburg/ESO surveys. These same surveys have tended to target the 
diffuse Galactic halo, but dwarf galaxies, the inner halo, the Galactic Centre, or another site, may prove to be a more 
prosperous hunting ground in future years (see Sec. IV. A I and Sec. IVIip . 

In the closing statements of many reviews in astronomy, there are the obligatory requests for more theoretical 
and numerical work in support of a rich harvest of observations. But in the field of pre-galactic enrichment, the 
observations appear more limited, in part because we have not yet learnt to properly interpret the data. There is a 
pressing need for more computational work on many fronts. These include: 



• better stellar atmospheric models to improve abundance calibrations for metal poor stars; 

• a new generation of stellar evolution models (especially giants, sub-giants) with updated nuclear reaction net- 
works, diffusion and mixing; 

• self-consistent stellar models that take into account the possible binary nature of the first stars; 

• self-consistent supernova models that properly treat the explosion mechanism, fallback, rotation and mixing; 

• better fine-scale and coarse-scale mixing theories and algorithms for metal transport in turbulent media; 

• binary star mass transfer (zero metallicity) with a proper treatment of Roche-lobe overflow, slow and fast winds; 

• higher resolution hydrodynamical models of the formation of the first stars and star clusters. 

Just where to look for the most ancient and/or second generation stars (and possibly metal free stars) is a hot topic 
at the present time with some groups favouring t he inner parts of ga.laxies (e.g. Tumlinson, 2010; White fc S pringel 
(2000)), and others finding a wide spread in radius ( Brook et al . 2007 : Salvadori et al . 2010t Scannapieco et aZ.l l2006). 
More advanced numerical simulations may well provide better guidelines on where to look and even the specifics of the 
orbit parameters of the oldest stars (e.g. Tumlinson, 2010). In this respect, we look forward to the Gaia astrometric 
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satellite which is due to launch in the next few years. This will provide exceptional kinematic information for up to a 
billion stars in the Galaxy. We recommend that future simulations attempt to predict the energy- angular momentum 
space (e.g. Helmi, 2008) expected for these ancient stars. 

Regardless of the simulations , it is important to keep in mind that metal-rich gas is observed in quasar/AGN 
spectra to the hig hest redshifts (iHam ann and Ferland|j I1999I) consistent w ith detections of dust and molecular gas at 
these early times ( Cox et al. . 2002; Klamcr et at, 20041: lYun et al. . 2000f ). This reflects the fact that the dynamical 
timescales in the cores of galaxies are very short lO^yr). Thus ancient stars with solar abundances may well exist 
in the cores of galaxies, assuming these have not been flung out into the disk by the formation of a central bar (e.g. 
Minchev and Famaey, 2009; see Sec. IV.B|) . 

Precisely when the first star clusters formed is unclear. Interestingly, the difficulty in iden tifying metal-free s t ars at 



the present epoch may already provide a constraint on the primordial initial mass function ( Scannapieco et al. . 20061 ) 



and the possible existence of star clusters. As shown in Table IIIIl a star must be less massive than the Sun to have 
survived a Hubble time. If the star formed at a redshift as late as 2 ^ 1, its mass must be less than 0.9 — 1.1 M-q, 
depending on the metallicity, in order to have reached the giant branch today. If these stars did not exist at those early 
epochs due to a truncated IMF, then we must infer the progenitor yields from later generations of stars. Interestingly, 
because 'clustered' abundance signals will be unaffected by mass transfer in binaries, these may give us a truer reading 
of the IGM metallicity at early times. 

A now approach to exploring the evolution of the early initial cluster mass function has bee n proposed by 



Bland-Hawthorn et aLl (l2010ah . Stars born in clusters are extremely uniform in their chemical elements ( Bland- Hawthorn et al 
2010bl : iDe Silva et all l2007bh . This should lead to clumping in the abundance plane, e.g. [Fe/H] vs. [a/Fe]. While 



the effect can be detected in 8m class data, a much richer return is expected in an era of extremely large telescopes. 
To this end, it will be necessary to equip these telescopes with wide-field multi-object spectrographs that operate at 
high spectroscopic resolution (i? >20 000). Such an instrument can be profitably targetted at dwarf galaxies, the 
Galactic bulge or populations identified by the Gaia satellite. 
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